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The  use  of  new  advanced  solid  rocket  propellant  constituents,  such  as,  e.g., 
boron  as  HEDM  and  ammonium  nitrate  as  replacement  for  halogen-containing 
oxidants,  requires  knowledge  of  their  combustion  characteristics.  To  this  end 
experimental  rate  coefficient  k  measurements  on  individual  important  A1  and  B 
species  reactions  were  made  in  the  300-1800  K  temperature  range  using  a  hig  - 
temperature  fasi-flow  reactor.  Temperature  dependent  rate  coefficmnt  expressions 
kiTl  =  AT(super  n)exp(-E/RT)  are  found  for  the  reactions  between  AlO  and  CH(sub  4), 
O^b  2)  A1  and  NXsub  2)0:  AlCl  and  N(sab  2)0.  CO(sab  2),  SO(sub  2);  BCl  and  Msub 
2)0  SO(sub  2);  BO  and  CO(sub  2).  To  allow  semi-empincal  predictions  for  further 
metal  oxidation  reactions,  two  semi-empirical  theories  have  been  developed  and 
utilized.  These  are  Semi-Empirical  Configuration  Interaction  theory  for  prediction  of  E 
from  one  species  to  another  and  Metals-Transition  State  Theory  for  prediction  of  k(T) 
from  k  at  one  temperature.  Dual  military-civilian  applications  of  this  work  are 
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I.  RESEARCH  OBJECTIVES 


Advances  in  solid-fuel  propellant  technology  for  the  late  nineteen  nineties 
depend  critically  on  an  understanding  and  quantitative  description  of  a  number  of 
combustion  processes.  These  include  the  combustion  of  metallic,  particularly  boron, 
HEDMs  and  environmentally  clean  A1  and  B  propellants,  containing,  e.g.  ammonium 
nitrate  instead  of  perchlorate.  In  this  context  the  present  work  had  several 
objectives.  Unique  high-temperature,  fast-flow  reactors  HTFFR,  which  allow 
measurements  of  the  kinetics  of  isolated  reactions  in  the  300-1900  K  temperature 
range,  were  employed  to  provide  a  kinetic  data  base  for  metallic  species  reactions  in 
propulsion  system  models.  The  results  obtained  enhance  the  understanding  of  the 
chemistry  involved,  which  in  turn  allowed  making  predictions  for  further  reactions. 
This  enables  faster  model  development  than  experimental  work  by  itself  can  provide. 
To  this  end  establishing  correlations  in  series  of  similar  reactions  and  understanding 
the  fundamental  basis  for  these  was  emphasized.  Where  appropriate  dual  military- 
civilian  use  aspects  of  the  work  have  been  emphasized  (see  especially  Section  III  A). 


II.  RESULTS 


A.  Cumulative  Chronological  List  of  Publications 
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1 .  Arthur  Fontijn,  Ed.  Gas-Phase  Metal  Reactions,  North-Holland,  Amsterdam, 
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2.  David  P.  Belyung,  Arthur  Fontijn  and  Paul  Marshall,  "Gas-Phase  Reactions 
Between  Hydrocarbons  and  Metal  Oxides.  The  AIO-1-CH4  reaction  from  590 
to  1380  K",  The  Journal  of  Physical  Chemistry,  97.,  3456-3459  (1993). 

APPENDIX  C 


3.  Peter  M.  Futerko  and  Arthur  Fontijn^  "Experimental  and  Transition- 

State  Theory  Studies  of  the  Gas-Phase  Reactions  of  AlCl  with  N2O,  CO2, 
and  SO2,"  The  Journal  of  Physical  Chemistry.  97.  7222-7227  (1993). 
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(1993). 
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Barriers  for  Series  of  Exothermic  Homologous  Reactions  IV.  Comparison  of 
Measurements  to  Theory  for  Reactions  of  s^p^  Atoms  with  N2O,  The 
Journal  of  Chemical  Physics  XQTli  155-160  (1995). 
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Institute,  1993,  University  Microfilm  No.  9325242. 


ni.  SIGNIFICANCE  OF  THE  RESULTS 

Here  we  discuss  the  significance  of  the  papers  given  in  the  preceding  section  in 
the  chronological  sequence  given  there. 

A.  General  and  Military-Civilian  Dual-Use  (Publications  1  and  2) 

The  general  significance  of  Publication  1,  should  be  clear  from  its  Preface  given. 

It  has  become  very  widely  quoted. 

A  few  dual  military-civilian  aspects  are  illustrated  here.  Publication  2  is  one 
example  of  this;  another  is  the  use  of  our  unique  family  of  reactors  developed  for 
wide  temperature  range  studies  on  isolated  reactions.  The  development  of  these 
reactors  was  prompted  by  the  needs  of  the  military  for  kinetic  data  on  metal 
combustion.  With  support  of  NSF  we  are  now  also  applying  the  reactors  to  waste 
incineration  problems.  In  that  area  we  showed  the  presence  of  several  paths  for 
chromium  oxidation.  1  As  chromium  is  in  its  highest  oxidation  state  highly  toxic,  but 
relatively  harmless  in  lower  oxidation  states,  knowledge  of  the  various  reaction 
pathways  is  needed  to  hinder  the  complete  oxidation.  An  understanding  of  high- 
temperature  reactions  of  metals  is  important  to  many  more  applications.  These 
include  decreasing  slag  formation  in  coal  combustion,  preventing  hot  corrosmn  in  jet 
engine  turbines,  alleviation  of  soot  formation  in  all  sorts  of  combustion  devices,  as 
well  as  new  chemical  synthesis  processes  for  refractory  materials. 

The  direct  conversion  of  methane,  i.e.  natural  gas,  to  liquid  fuels  and  chemicals 
of  commercial  interest  is  a  long-sought  goal.  A  major  approach  being  followed  is  to 
pass  methane  over  solid-state  catalysts  consisting  of  mixtures  of  metal  oxides  with 
dopants.  The  first  step  in  this  process  is  to  abstract  a  hydrogen  atom  from  methane 
CH4  to  produce  a  highly  reactive  methyl  radical  CH3  to  form  C-C  bonds.  To 
understand  these  surface  reactions  better,  gas-phase  studies  of  the  reactions  between 
methane  and  metal  oxides  are  desired.  Theoretical  studies  have  suggested  that  such 
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abstraction  reactions  should  occur  at  elevated  temperatures.  Because  of  our 
experience  with  AlO  from  studies  for  rocket  propellants,  we  decided  to  use  that  oxide 
to  experimentally  investigate  this.  Publication  2  provides  the  first  experimental 
confirmation,  but  also  shows  that  at  least  in  this  instance  the  theoretically  predicted 
mechanism  is  not  dominant.  In  general,  such  gas-phase  studies  provide  essential 
input  for  gas-solid  catalytic  reactions. 

B.  k(T)  Predictions  for  Model  Development  (Publications  3-5) 

Combustion  models,  whether  for  the  plume  or  chamber,  require  reliable  kinetic 
input  data  at  realistic  temperatures.  The  data  needed  are  rate  coefficients  k(T)  and 
mechanisms.  Measurements  with  our  HTFFR  (high-temperature  fast-flow  reactor) 
and  MHTP  (metals  high-temperature  photochemistry)  techniques  are  providing  such 
information  for  metallic  species.  However,  experiments  alone  are  a  rather  slow  way 
to  obtain  information  on  all  the  reactions  needed.  It  is  therefore  desirable  to  develop 
predictive  techniques.  Fully  theoretical  a  priori  techniques  have  not  yet  reached  the 
state  of  making  meaningful  predictions.  We  have  therefore  chosen  to  develop  semi- 
empirical  techniques  for  metallic  species.  These  predictions  allow  either  extrapolation 
of  our  measured  data  from  one  species  to  another  (SECI  theory),  or  for  one  reaction, 
from  one  temperature  to  a  wide  range  of  temperatures  (MTST  theory).  Another  way 
of  expressing  the  scope  of  these  theories  is  that  in  rate  coefficient  expressions 

k(T)  =  AT"  exp(-E/RT)  (1) 

SECI  predicts  the  exponential,  MTST  the  pre-exponential  part.  The  theories  are  semi- 
empirical  in  that  for  SECI  the  E  factor  (activation  barrier)  of  a  similar  reaction  needs 
to  be  known  and  that  for  MTST  an  experimental  rate  coefficient  measurement  at  one 
temperature  is  required. 

The  left  hand  side  of  Figure  1  illustrates  SECI  the  right  hand  side  MTST,  each  of 
which  will  now  be  discussed. 

1 .  Semi-empirical  Configuration  Interaction  (SECI)  Theory 

The  E  factor  in  eq.  (1)  is  the  difference  in  the  potential  energy  of  the  activated 
complex  (transition  state),  formed  from  the  reactants,  and  the  sum  of  the  potential 
energies  of  the  reactants.  It  represents  the  highest  potential  energy  on  the  path  from 
reactants  to  products.  It  can  be  calculated,  using  resonance  theory,  from  the 
properties  of  three  hypothetical  resonating  activated  complexes,  i.e.  that  resulting 
from  purely  ground  state  interaction  for  which  the  barrier  is  q,  the  interaction 
involving  the  lowest  excited  state  of  the  atomic  reaction  partner,  and  ionic  interaction. 
The  wave  function  of  the  actual  complex  resulting  from  these  structures  is  given  by 

T  =  ci'Fi  +  C2T2  +  C3T3  (2) 

where  the  c  factors  give  the  contributions  of  each  structure.  The  actual  T  will  be  that 

for  which  E  =  J  'F*HTdx  is  at  a  minimum.  The  condition  can  be  expressed  as  a  set  of 
simultaneous  linear  equations,  ultimately  resulting  in 
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(3) 


aq  +  bq^  +  cqE  +  dq^E  +  fE  +  gE^  +  hE^  -  0 

where  a,  b,  c,  d,  e,  f,  g,  h  are  constants  derived  from  tabulated  physical  properties  of 
the  individual  reaction  partners. 

E  and  q  can  not  simultaneously  be  calculated  from  first  principles.  For  a  series 
of  similar  reactions,  e.g.  the  boron  group  atoms  with  N2O  (see  the  figure)  we  therefore 
obtain  E  for  one  atomic  partner.  This  is  done  by  fixing  n  for  the  series,  which  allows 
finding  E  from  eq.  (1)  used  as  the  fitting  expression  for  the  experimental  k(T) 
measurements.  Eq.  (3)  then  yields  q,  which  according  to  the  Hirschfelder  theory  is 
the  same  for  all  members  of  the  series.  E  for  the  other  members  of  the  series  can 
now  be  calculated  from  eq.  3  using  q  and  the  constants  appropriate  for  each  of  the 
atoms. 

The  SECI  figure  shows  E  plotted  against  I.P.  (the  ionization  potential  related  to 
the  3rd  structure)  plus  P.E.  (the  promotion,  or  excitation,  energy  representative  of  the 
2nd  structure).  It  shows  good  correlation  with  IP  +  PE  and  excellent  agreement 
between  theory  and  experiment  obtained  for  three  groups  of  atom  reactions  with  N2O. 
The  left  hand  line  is  for  the  one-electron  alkali  metal  group  and  is  of  interest  for 
example  to  hot  turbine  combustion  corrosion  problems.  The  middle  line  is  for  two 
electron  atoms,  and  contains  metals  of  interest,  e.g.  to  toxic  emission  problems.  The 
right  hand  line,  for  two  s  one  p  electron  atoms,  has  most  recently  been  obtained  and 
includes  metals  of  most  importance  to  rocket  propulsion. 

2 .  Metals  Transition  State  Theory  (MTST) 

As  shown  in  Fig.  1,  by  taking  one  experimental  point  (the  solid  symbol)  at  mid¬ 
temperature  range  for  the  reaction  of  interest,  we  can  predict  the  shape  of  the 
Arrhenius  plot  in  good  agreement  with  measurements,  thus  reducing  the  work 
required  for  data  taking. 

In  general  transition  state  theory  (TST)  k(T)  is  given  by: 


k(T)  = 


keT  q^ 
h  qABqCD 


exp(-AE’‘/RT) 


(4) 


where  ke  and  h  are,  respectively,  the  Boltzmann  and  Planck  constants,  and  qAB  and 
qCD  are  the  partition  functions  of  the  reactants,  q’*  and  AE^are,  respectively,  the 
partition  function  of  the  transition  state  and  its  energy  above  the  reactant  ground 
states,  both  of  which  are  unknown.  By  using  simple  ab  initio  methods  for  estimating 
and  equating  eq.  (4)  with  an  experimental  rate  coefficient  at  one  temperature,  AE* 
can  be  found. 


The  basic  TST  approach  was  developed  primarily  for  reactions  of  more 
commonly  studied  reactants,  such  as  hydrocarbons.  There,  the  transition  state 
structure  and  hence  q’*  can  be  found  readily  by  comparison  to  known  stable 
compounds,  which  information  is  generally  lacking  for  refractory  species.  The 
advance  is  that  we  have  made  the  method  useful  for  such  species  as  well. 
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C  AIO  and  BO  Oxidation  Reactions  (Publication  6) 

Oxidation  reactions  of  A1  and  B  species  are  essential  input  to  various  rocket 
combustion  models.  In  further  experimental  work,  we  have  thoroughly  studied  the 
AIO  oxidation  by  O2  (Publication  6)  and  obtained  a  significant  preliminary  result  for 
BO  oxidation  (not  published). 

The  AIO  +  O2  reaction  system  was  originally  thought  to  proceed  as  a  simple 
abstraction  reaction  leading  to  AIO2  +  O.^  However,  in  recent  years  we  have  found 
several  examples  of  reactions,  considered  in  the  literature  as  simple  abstraction  to 
proceed  by  addition  or  insertion  at  lower  temperatures,  with  abstraction  becoming 
more  significant,  but  not  necessarily  dominant,  at  higher  temperatures.  Examples  of 
this  are  publication  2  and  Ref.  1.  These  observations  and  the  vast  improvements  in 
HTFFR  technique  since  then,^  prompted  this  new  investigation.  For  temperatures 
below  about  1000  K  the  AIO  consumption  now  is  found  to  be  pressure-dependent, 
with  a  negative  temperature  dependence;  log  k  (305-1010  K)  =  -25.36-1.69  log  (T/K) 
cm6molecule-2s-i.  At  higher  temperatures  the  system  becomes  pressure- 
independent  and  the  AIO  consumption  rate  coefficients  increase,  k(1200-1690  K)  — 
7.7xl0-i0exp  (-10008  K/T)  cm^molecule-is-^. 

The  following  mechanism  describes  these  findings: 


AIO  +  O2 


AlO^  +0 


AIO3 


(5) 


Here,  the  dotted  lines  indicate  the  paths  unlikely  to  contribute  significantly.  Thus,  the 
reaction  is  thought  to  proceed  through  an  intermediate  AIO3  complex,  which  can  be 
collision-stabilized  or  dissociate  to  AIO2  +  O.  The  latter  process  is  96  +.  40  kJ  mol 
endothermic. 


It  is  interesting  to  observe  that  in  a  room  temperature  study  of  the  AIO  +  CO2 
reaction,  competition  between  an  addition  and  abstraction  reaction  was  observed.3 
This  suggests  that  in  that  system  abstraction  will  become  dominant  at  lower 
temperatures  than  for  AIO  +  O2.  The  O2  reaction  should  be  primarily  of  interest  for 
plume  models,  while  the  CO2  reaction  is  of  both  plume  and  chamber  interest. 

The  AIO  +  O2  behavior  contrasts  with  that  of  the  iso-electronic  exothermic 
BO  +  O2  reaction  system,  for  which  no  pressure  dependence  was  found  over  the  300- 
1000  K  range,  but  again  a  negative  temperature  dependence  was  observed;^ 
k(T)  =7.0xl0-l2exp  (259  K/T)  cm3molecule-is-i.  This  suggests  a 


bo  +  02;!bo3  -^bo2+o  (6) 

mechanism.  We  have  made  a  series  of  BO  oxidation  measurements  in  the  HTFFR 
apparatus.  In  this  environment  accurate  rate  coefficients  could  not  be  obtained, 
apparently,  due  to  BO  wall-dimerization.  We  are  about  to  continue  this  work  in  one  of 
our  MHTP  (metals  high-temperature  photochemistry)  reactors,  where  such  problems 
should  not  arise. ^  Nonetheless,  the  HTFFR  work  has  already  yielded  some  interesting 
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information.  Our  results  on  BO  +  O2  between  400  and  1000  K  are  in  essential 
agreement  with  the  earlier  measurements  from  an  MHTP-like  reactor.^  For  the 
hypothesized  reaction  BO  +  CO2  BO2  +  CO  a  slow  pressure-dependent  process  was 
observed  at  room  temperature,  but  no  measurable  BO  consumption  occurred  from  780 
to  1200  K.  This  indicates  rate  coefficients  less  than  about  lxlO-^5  cm^molecule'^s'^. 
This  is  four  orders  of  magnitude  lower  than  found  for  the  BO  -1-  O2  reaction.  It  further 
is  two  to  three  orders  of  magnitude  lower  than  the  guestimated  values  in  the  current 
models  for  boron-hydrocarbon  slurry  (air-breathing)  combustion,  where  it  is  indicated 
to  be  a  highly  sensitive  reaction.  Those  models  should  now  be  reconsidered. 


IV.  PROFESSIONM.  PERSONNEL 

For  two  Ph.D.  students,  P.M.  Futerko  and  D.P.  Belyung  this  grant  and  the  parallel 
AASERT  grant  were  their  principal  projects.  Additionally,  Dr.  P.N.  Bajaj  and  Mr.  W.F. 
Flaherty  helped  with  the  work,  as  did  two  M.Sc.  students  M.G.  Santana  and  G.T. 
Dalakos,  and  an  undergraduate  A.M.  Zinacola.  My  former  students.  Professor  P. 
Marshall  and  Dr.  A.G.  Slavejkov  contributed  significantly  to  publications  2  and  4, 
respectively. 

V,  CUMULATIVE  CHRONOLOGICAL  LIST  OF 
PRESENTATIONS  AND  OTHER  INTERACTIONS 

We  presented  papers  and  seminars  in  which  results  of  our  AFOSR-sponsored 
work  were  discussed,  at  the: 

1.  HEDM  Contractors  Conference,  Lancaster,  CA  (April  1992). 

2.  Russian  Academy  of  Sciences,  Chernogolovka,  Russia  (May  1992). 

3.  American  Chemical  Society,  National  Meeting,  Washington,  DC  (August  1992). 

4.  Department  of  Chemistry,  University  of  Auckland,  New  Zealand  (March  1993). 
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PREFACE 

The  study  of  reactions  of  metallic  species  in  the  gas  phase  impinges  on 
a  variety  of  basic,  technical  and  socio-technical  fields.  Spread  over  the 
Lr^ntlv  developed  world  there  are  a  sizeable  number  of  scientists  and 
TnSZk  occupied  with  gas-phase  metal  reactions.  \et.  there  has  been 
Sue  interaction  among  them;  apparently  no  book  or  major  meeting  has 
SevSuslv  been  dedicated  to  the  subject.  A  number  of  conversations  I  had 
with  colleagues  on  both  sides  of  the  Atlantic  dunng  a  recent  sabbatical 
rlearlv  brought  out  the  need  to  try  to  bundle  our  knowledge.  The  present 
volume  has  come  about  as  a  response  to  this^  It  bnnp  together  overviews 
frormost  of  the  groups  concerned  with  fundamental  research  both  of  an 
Sn^nmental  and  theoretical  nature.  Still,  much  of  this  work  has  been 
"tubulated  bv  very  practical  problems.  The  individual  authors  show  these 
ronnections  'To  further  emphasize  this  aspect,  a  few  of  the  final  chapters 
pnmanly  address  such  problems,  but  show  what  the  research  needs  are.  or 
how  fundamental  information  can  be  used  toward  their  :5olution.  The 
introductory  chapter  outlines  the  connection  between  the  various 
contributions  and  mentions  significant  turther  efforts  and  applications. 

If  one  compares  the  present  field  to  other  areas  of  chemical  reaction 
-tud'C'  such  as  hydrocarbon  oxidation,  it  is  striking  how  little  theoretical 
-iTort  there  has  been.  In  fact,  several  of  the  theory  chapters  in  this  volume 
came  about  pnmanlv  because  the  plans  for  it  stimulated  the  research. 
.Another  gap  that  shows  clearly  is  that  between  the  kinetic  and  dynamic 
aspects;  much  work  needs  to  be  done  before  a  proper  synthesis  of  the  results 
obtained  can  be  achieved.  Many  of  the  chapters  indicate  current 
shortcomings  and  potential  areas  where  important  contnbutions  can  be 
made  This  book  therefore  not  only  summarizes  most  of  the  present 
knowledge,  but  also  provides  a  proper  starting  point  for  far  more  extensive 
further  research  and  interdisciplinary  collaboration. 

The  birth  of  this  volume  has  proceeded  in  stages.  The  authors 
submitted  preliminary  chapters,  which  were  then  presented 
at  a  symposium  held  before  the  Physical  Chemistry  Division  of  the  Amenta 
Chemical  Society,  m  New  York  in  August  1991.  Following  meeting 
many  of  the  chapters  we.-e  revised  and  all  were  '^en  refereed  and  usually 
further  improved.  Thus  this  work  already  includes  the  benefit  of  much 
interaction  on  an  international  scale. 
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The  «tudv  of  reactions  of  metallic  species  in  the  gas  phase  impinges  on 

of  basic  technical  and  socio-technical  fields.  Spread  over  the 

ru^r^ntW  developed  world  "  there  are  a  sizeable  number  of  scientists  and 
pnmneers  occupied  with  gas-phase  metal  reactions.  \et.  there  has  been 
litUe  interaction  among  them;  apparently  no  book  or  major  meeting  has 
nievuoSSv  been  dedicated  to  the  subject.  A  number  of  conversations  I  had 

^vafh  colleagues  on  both  sides  of  the  Atlantic  dunng  a  recent  sabbatical 
Aitn  coiieag  knowledge.  The  present 

Sume  h^?le  ab^^^^^^  response  to  this^  It  brings  together  overviews 
colume  groups  concerned  with  fundamental  research,  both  of  an 

mental  and  theo^redca  nature.  Still,  much  of  this  work  has  been 
expenmental  and  tneore^i^^^  problems.  The  individual  authors  show  these 

not'  To'fur^her  this  aspect,  a  few  of  the  final  chapters 

pnmanlv  address  such  problems,  but  show  what  the  research  needs  are,  or 
Lw  fundamental  information  can  be  used  toward  their  =olutiom  T 
ntrodlictorv  chapter  outlines  the  connection  betw^een  the  various 
contnbutions  and  mentions  significant  further  efforts  and  applications. 

If  one  compares  the  present  field  to  other  areas  of  chemical  reaction 
>;tudie-  such  as  hydrocarbon  oxidation,  it  is  striking  how  little  theoretica 
Srt  the%  has  been.  In  fact,  several  of  the  theory  chapters  in  this  volume 
came  about  pnmanlv  because  the  plans  for  it  stimulated  the  research, 
Cther  gap  that  shows  clearly  is  that  between  the  kinetic  and  dynamic 
aspects-  much  work  needs  to  be  done  before  a  proper  synthesis  of  the  results 
obtained  can  be  achieved.  Many  of  the  chapters  indicate  current 
shortcomings  and  potential  areas  where  important  contnbutions  can  b 
made  Thi'  book  therefore  not  only  summarizes  most  of  the  present 
r„„  .'l.dge:b„ T.lso  prev.dbs  a  prop*.-  atartms  pent  for  far  mora  ea.ens.ve 
further  research  and  interdisciplinary  collaboiation. 

The  birth  of  this  volume  has  proceeded  in  stages.  The  authors  ^  ^ 
submitted  preliminarv  chapters,  which  were  then  presented  and  discussed 
at  a  symposium  held  before  the  Physical  Chemistry  Division  of  the  .Amencan 
Chem'icafsocietv  in  New  York  in  August  1991.  Following  the  meeting 
manTof  the  Chapters  w^re  revised  and  all  were  then  refereed  and  usually 
further  improved.  Thus  this  work  already  includes  the  benefit  of  much 
interaction  on  an  international  ::cale. 
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A  recent  ab  inito  study  suggested  that  gaseous  metal  oxides  can  dirertly  abstrart  H  atoms  from  hydrocarbons, 
but  with  considerable  barriers.  To  try  to  confirm  the  occurrence  of  such  reactions,  the  title  reaction  has  been 
studied  in  a  high-temperature  fast-flow  reactor.  The  data  are  well  described  by  the  fitting  expression  fc(590- 
1380  K)  =  9.6  X  10-^’(r/K)^’‘  exp(2468  K/7)  cm^  molecule*'  s*',  with  2a  precision  limits  varying  with 
temperature  from  ±4%  to  ±  1 2%.  The  corresponding  2a  accuracy  limits  arc  about  ±25%.  Comparison  of  this 
result  to  further  ab  initio  and  conventional  transition-state  theory  calculations  suggests  that  dir^  abstraction 
indeed  can  occur  at  the  hi^Jer  temperatures  but  that  one  or  more  other  channels,  possibly  involving  AlO 
insertion  into  a  C-H  bond,  dominate  in  the  initial  attack  step.  The  potential  implications  of  this  work  for 
catalytic  conversion  of  methane  to  higher  hydrocarbons  are  considered. 


Introduction 

In  a  recent  ab  initio  study  Borve  and  Pettersson*  predicted 
that  some  gaseous  metal  oxides,  LiO,  MgO,  and  AlO,  can  directly 
absuact  an  H  atom  from  methane  to  produce  methyl  radicals. 
Such  exothcnnic  reactions  have  not  been  observed  in  the  gas 
phase,  and  their  studies  suggested  these  to  have  considerable 
barriers,  25-^7  kJ  mol*' .  However,  solid  metal  oxides  have  been 
found  to  be  effective  catalysts  for  H  abstraction  from  methane, 
leading  to  C2H4  and  C2H6  in  significant  yields,^  a  process  of 
considerable  commercial  promise.  The  basis  for  this  catalytic 
reaction  has  been  strongly  suggested  to  be  the  surface  reaction 
of  the  O"  ion  to  remove  a  hydrogen  atom,  a  process  which  appears 
to  proceed  to  a  significant  degree  only  above  about  1000  K.^ 
Theoretical  studies  of  cluster  models  of  the  M0O3  oxide  structure 
have  supported  this  hypothesis.^  If  the  occurrence  of  the  gas- 
phase  processes  could  be  experimentally  confirmed,  it  would  not 
only  establish  a  new  class  of  reactions  but  also  could  lead  to  a 
better  understanding  of  the  heterogeneous  catalytic  processes. 
We  therefore  decided  to  investigate  one  of  these  suggested 
reaaions: 

AlO  +  CH4  ^  AlOH  +  CH3  ( 1 ) 

for  which  A/T*  =  -27  kJ  mol-'.^*^  While  this  reaction  has  a 
higher  predicted  barrier  than  those  of  the  other  oxides,'  we  have 
had  extensive  experience  with  the  study  of  the  gas-phase  reactions 
of  AlO  at  temperatures  up  to  1600  K  in  HTFFRs  (high- 
temperature  fast- flow  reactors)^  and  therefore  selected  this  oxide. 

Techniqiie 

The  HTFFR  facility  used  in  this  work  and  the  data  handling 
procedures  have  been  previously  described. Briefly,  a  vertical 
mullitc  (McDanel  MV-30)  reaction  tube  (60  cm  long,  2.2  cm 
i.d.)  is  heated  to  the  desired  temperature  by  columns  of  resistively 
heated  SiC  rods  inside  an  insulated,  water-cooled  vacuum  housing. 
A1  vapor  is  produced  by  resistively  heating  an  Al-wcttcd  W  coil 
in  a  flow  of  Ar  carrier  gas.  A  trace  of  oxidizer,  0.5%  N2O  or  5% 
O2  in  Ar,  at  ^=4  X  10“^%  of  the  main  Ar  flow,  is  passed  through 
a  Pt  side  tube  to  a  location  just  downstream  from  the  vapor  source 
to  rapidly  convert  A1  to  AlO.*-'^  Further  downstream  CH4/ Ar 
mixtures  at  flow  rates  of  4-8%  of  the  main  Ar  flow  rate  are 
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introduced  through  a  movable  inlet.  Rate  coefficient  measure¬ 
ments  are  made  under  pseudo-first-order  conditions,  [AlO]  « 
[CH4],  in  the  stationary  inlet  mode,"  at  reaction  zone  lengths 
of  10  or  20  cm. 

Relative  AlO  concentrations  were  monitored  by  laser-induced 
fluorescence  (LIF)  using  a  pulsed  Lambda  Physik  EMG  101 
excimer/FL  2002  dye  laser  in  combination  with  a  KDP  doubling 
crystal.  Two  transitions  B^2^X^S  and  were  used.  In 

the  former,  AlO  is  pumped  on  the  464.8-nm  (1 ,0)  band  and  the 
fluorescence  is  observ^  through  a  482-nm  (20-nm  fwhm) 
interference  filter,  i.c.,  mainly  at  the  486.6-nm  (1,1)  band. '  ^  The 
B^r-X^S  system  could  only  be  used  up  to  1200  K  because  of 
interference  by  the  background  radiation  from  the  reactor  walls. 
For  measurements  above  1200  K  and  some  checks  at  lower 
temperatures,  the  C^S-X^S  (0,0)  transition  at  302.2  nm'^  was 
used  in  combination  with  a  301-nm  (1 1-nm  fwhm)  interference 
filter.  The  fluorescence  was  detected  by  an  EMI  9813QA 
photomultiplier  tube,  connected  to  a  Data  Precision  Analogic 
6000/620  100-MHz  transient  digitizer. 

The  gases  used  were  99.998%  Ar  from  the  liquid  Ar  (Linde), 
0.5%  N2O  (99.99%)  in  Ar  (99.999%)  from  Matheson,  5%  O2 
(99.99%)  in  Ar  (99.995%)  from  Scott,  and  CH4  (99.99%)  from 
Matheson, 

Results 

Plots  of  In  [AlO]  relative  vcrsus  [CH4]  yield  straight  lines  with 
slopes  -kt,  where  /  is  reaction  time.  For  each  individual 
measurement,  k  and  Ok  were  determined  by  a  weighted  linear 
regression."  '^  The  k  values  and  the  experimental  conditions 
under  which  they  were  obtained  are  summarized  in  Table  I.  The 
measurements  may  be  seen  to  be  independent  of  the  following: 
pressure  varied  from  13.1  to  63,3  mbar,  corresponding  to  toul 
concentrations  [M]  from  8,0  x  lO'^*  to  6.4  X  10'“'  molecules 
cm'^;  average  gas  velocity  v  varied  from  1 2  to  98  m  S"';  reaction 
zone  length  selected  at  10  or  20  cm;  initial  fluorescence  intensity 
F  (a  measure  of  initial  [AlO])  varied  from  9  to  65  in  arbitrary 
units.  The  temperature  range  covered  is  590-1 380  K.  The  lower 
limit  was  determined  by  the  heating  effect  of  the  A1  vaporizer, 
and  the  upper  limit  by  the  dissociation  of  CH4.  An  Arrhenius 
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TABLE  I;  Summary  of  Rate  Coefficient  Measurements  on  the  AlO  +  CH.  Reaction*  - 


reaction  zone 
length  (cm) 


/>(mbar)  [M]  (10*'  cm  -^)  [CH4]mu»  (10*'  cm  ■) 


t^(ms 


r(K)  it  ±  (cm^  molecule  *  s  ') 


20 
20 
20 
20 
20 
20 
20 
20 
20 
20 
20 
20 
20 
20 
20 
20 
20 
10 
10 
10 
10 
20 
20 
20 
20 
10 
10 
10 
10 
20 
20 
20 
20 
20 
10 
10 
10 
10 
20 
20 
10 
10 
10 
10 
10 
20 
20 
20 
20 
20 
20 
20 
20 
20 
20 
20 
20 
20 
20 
20 
10 
20 
20 
20 
20 
20 
20 
20 
20 
10 
10 
20 
20 


26.7 

26.5 

43.3 
20,0 
20.0 

15.5 

15.5 

13.1 

17.2 

17.2 
28.9 

19.2 

16.7 
20.1 

36.7 

30.4 

32.8 

32.8 

36.3 

19.6 
19.6 
19.6 

19.6 

31.7 

33.6 

34.3 
34.3 
34.3 

34.3 

34.4 

34.4 

24.1 

24.9 
24.9  . 

25.1 

25.7 
36.0 

35.7 
16.0 
16.0 
16.3 

16.3 

21.2 

39.5 

39.5 

34.1 
34.0 
34.0 

53.1 

24.5 

24.5 

14.3 

14.3 

63.3 

15.3 

15.3 

21.9 
56,0 
13.7 

13.7 

14.1 

44.5 

44.7 

41.7 

41.6 

53.9 

54.3 
54.3 
17.5 

24.1 
24.1 
24.0 
24.0 


2.1 

2.1 

3.3 
1.5 

1.5 
1.0 
1.0 
0,8 
1.1 
1.1 
1.9 
1.2 
1.1 

1.6 
2.6 
2.1 
2.2 

2.3 
2.5 

1.4 

1.4 
1.3 

1.3 

2.3 

2.5 

2.6 
2.6 
2.6 
2.6 
2.5 
2.5 
2.1 
2.5 
2.5 
2.5 
2.5 

3.3 

3.3 
1.2 
1.2 

1.3 

1.3 

2.4 
4.3 

4.3 
4.1 
4.0 

3.9 

5.3 

2.9 
2.9 
1.7 

1.7 

6.4 

1.4 

1.4 

1.8 
4.8 
1.2 
1.2 
1.2 

2.3 

2.5 

4.5 

4.6 

5.3 

5.6 

5.7 
l.O 

1.4 
1.3 
1.3 
1.3 


9.4 

9.4 
10.3 

4.8 

4.8 
5.0 
5.0 

4.4 

3.8 
3.8 

4.7 

5.8 

4.9 

9.9 

4.4 
5.8 

5.1 

8.5 

9.4 

6.4 

7.7 
5.0 
5.0 
7.0 

7.5 

9.6 
9.5 
7.4 

7.4 

5.7 

5.7 

13.7 

18.4 

18.1 
18.6 

18.8 

15.5 
15.2 

5.6 
4.5 

6.9 
8.1 

14.0 

17.1 

17.2 
21.1 

17.1 

16.9 

6.2 

20.6 
20.0 
11.0 

11.1 

9.9 

8.7 

8.8 

5.9 
4.7 
4.7 

4.7 
7.0 

1.1 

1.2 

7.4 

7.5 

6.8 

8.3 

8.5 

1.5 
2.0 
2.0 

1.3 
1.3 


17 

17 

24 

33 
24 
22 
20 
30 

47 
32 
38 
41 
64 

53 
30 
41 

38 
35 
32 
64 
63 

54 
53 

48 

34 

37 

39 
29 

24 

25 

23 

26 
28 

24 
28 
28 

52 

38 

53 

35 

39 
65 
34 


32 

33 
29 
62 
62 
98 
98 
90 
68 
68 
48 
56 
67 
46 

29 

33 
31 

30 
28 
51 

51 

52 
52 

37 
35 

34 
34 
44 
44 
46 
46 

38 
29 
29 
28 
28 

34 

35 
58 
58 
57 
57 
28 


913 

912 

928 

932 

935 

1083 

1083 

1135 

1098 

1102 

1103 

1139 

1144 

882 

1001 

1051 

1053 

1031 

1034 

1034 

1037 

1061 

1063 

985 

982 

961 
963 

962 
961 

1000 

1003 

827 

721 

734 

717 

734 

781 

795 

945 

953 

938 

942 

651 


9 

15 

660 

9 

15 

663 

42 

19 

602 

30 

19 

618 

25 

19 

628 

9 

21 

727 

15 

19 

603 

27 

20 

621 

33 

24 

598 

20 

23 

591 

11 

20 

718. 

34 

37 

780 

40 

37 

770 

25 

33 

893 

23 

20 

845 

20 

55 

846 

16 

56 

853 

24 

56 

866 

19 

30 

1374 

20 

28 

1380 

32 

13 

668 

45 

12 

650 

62 

14 

736 

65 

13 

707 

58 

13 

686 

40 

60 

1232 

37 

46 

1283 

37 

47 

1296 

38 

48 

1317 

39 

48 

1324 

4.42iC.33(-14)‘--^ 
4.93  ±0.41  (-14)-^ 
3.19  ±0.31  (-14)‘^ 

5.22  ±0.63  (-14)*^ 
6.78  ±0.93  (-14)^ 
1.67  ±0.21  (-13)^^ 
1.83  ±0.21  (-13)^ 
1.76  ±0.23  (-13)^ 
1.72  ±0.21  (-13)^ 
1.75  ±0.17  (-13)'^ 
1.09  ±0,12  (-13)'^ 
1.58  ±0.19  (-13)^ 

1.81  ±0.22  (-13)^ 
3,30  ±0.53  (-14)^ 

8.15  ±0.90  (-14)^ 
8.53  ±0.79  (-14)^ 

‘  1.11  ±0.08  (^13)^ 

8.82  ±0.97  (-14)^ 
8.62  ±0,81  H4)^ 
1.43  ±  0.16  (-13)^^ 

9.28  ±  1.07  (-14)-^ 

1.28  ±0.14  (-13)^ 

1.29  ±  0.15  (-13)'^ 
7.38  ±  0.53  (-14)^ 

7.22  ±0.55  (-14)^ 
7.91  ±0.69  (-14)^ 

8.18  ±0.66  (-14)^' 

7.11  ±0.63  (-14)-^ 

8.19  ±0.85  (-14)^ 
1.09  ±0.09  (-13)^ 
1.09  ±0.09  (^13)-^ 

2.10  ±0.21  (-14)^^ 

1.16  ±0,15  (-14)^ 
1,34  ±0.14  (-14)^ 
1.31  ±0.12  (-14)^ 
1.65  ±0.22  H4)^ 
2.67  ±0.32  (-14)^ 

3.12  ±0.39  (-14)^ 
7.05  ±  0.83  (-14)' 
7.84  ±0.97  (-14)' 
MO  ±0.12  (-13)' 
1.07  ±0.14  (-13)' 
1.44  ±0.18  (-14)' 
1.62  ±0.20  (-14)' 
1.99  ±0.14  (-14)' 
5.60  ±0.64  (-15)' 
9.95  ±  1.27  (-15)' 
1.02  ±0.13  (-14)' 

2.11  ±0.23  (-14)' 
7.78  ±0.91  (-15)' 

8.11  ±0.65  (-15)' 

6.71  ±0.77  (-15)' 
6.23  ±0.78  (-15)' 

1.72  ±0.20  (-14)' 
2,40  ±0.26  (-14)' 

2.43  ±0.31  (-14)' 
4.94  ±0.54  (-14)' 
3.02  ±0.31  (-14)' 
4.64  ±  0.53  (-14)' 

3.87  ±0.47  (-14)' 
5.81  ±0.75  (-14)' 
4.57  ±0.32  (-13)' 
3.49  ±0.34  (-13)' 
1.46  ±0.11  (-uy^ 
1.36  ±0.10  (-14/ 
1.71  ±0.11  (-14)/ 
1.68  ±0.14  (-14)/ 
1.45  ±0.11  (-14)/ 

3.44  ±0.37  (-13)' 
5.07  ±0.44  (-13)' 

5.87  ±0.50  (-13)' 
4.52  ±0.35  (-13)' 
5.21  ±0.44  (-13)' 


obtained  by  using  the  B-S-X-S  transition  and  a  quartz  ring  inlet. 
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plot  of  these  data.  Figure  1,  shows  strong  curvature.  Several 
different  expressions  were  tried  to  fit  the  In  k  versus  T~^  data. 
Using  Marquardt’s  method,’^  a  regression  fit  to  a  k(T)  =  A(T/ 
K)"  cxp(-BfT)  expression  yields 

A:(590-1380  K)  =  9.6  X  10-”(7'/K)’’‘  X 

cxp(2468K/7^  cm^  molecule"' s*'  (2) 

Equation  2  must  be  seen  as  a  fitting  expression  only,  as  the  numbers 
are  unrealistic  for  a  single  channel  reaction.  This  k{  T)  expression 
most  likely  indicates  competition  between  several  reaction 
channels.  We  next  attempted  a  similar  regression  for  a  double¬ 
exponential  expression  k{T)  =  A  expi-B/T)  +  C  expi-D/T), 
which  yielded 

Jt(590-1380  K)  =  2.0  X  lO"'"  exp(-8299  K/T)  + 

4.2  X  10"^^  cxp(-2402  K/T)  cm^  molecule'*  s'*  (3) 

This  fitting  expression  also  docs  not  seem  to  be  representative  of 
individual  actual  reaction  paths,  because  of  the  high  preexpo- 
ncntial  of  the  first  term.  This  term  appears  too  large  for  a 
metathesis  reaction  between  a  diatomic  radical  and  a  molecule.^ 
While  molecular  insertion  reactions,  e.g.,  of  singlet  CH2  and  SiH2, 
arc  known  with  A  factors  of  this  magnitude,  these  reactions 
proceed  with  an  activation  energy  close  to  zero.*^-’* 

Calculating  the  uncertainties  of  these  fitting  expressions  by 
combining*’  variances  and  covariances,  similarly  as  in  our  earlier 
studies,^  yields  precision  limits  of  ±4%  to  ±12%  for  cq  2  and  of 
±5%  to  ±12%  for  cq  3,  depending  upon  temperature.  Both 
expressions  lead  to  essentially  equal  accuracy  limits  of  ±23%  to 
±26%,  when  combined  with  a  ±10%  uncertainty  in  the  reactor 
flow  profilc“*2^  and  ±20%  for  potential  systematic  errors.  The 
simpler  expression  (2)  is  therefore  preferred  for  the  present 
temperature  range. 

Discussioo 

There  has  been  one  previous  attempt  to  observe  reaction  1 .  In 
that  room  temperature  study  no  reaction  could  be  detected,  and 
an  upper  limit  for  the  rate  coefficient  of  5  X  10'*^  cm^  molecule'* 
S'*  was  derived.*  Equations  2  and  3  suggest  much  lower  values, 
hence,  the  results  do  not  contradict  one  another. 

The  activation  energies,  i.c.,  the  local  slopes  of  the  Arrhenius 
plot.  Figure  1,  reach  the  67  kJ  moL*  calculated  by  Borve  and 
Pettersson*  only  at  the  highest  temperatures  investigated.  For 
example,  from  cq  2  these  energies  are  46, 59, 66,  and  71  kJ  mol  * 
at  1000,  1200,  1300,  and  1380  K,  respectively.  Also,  the  fitting 
expressions  suggest  that  the  reaction  follows  several  paths. 
Further  mechanistic  considerations  thus  appear  in  order.  To  this 
end  we  first  assess  the  likely  importance  of  H-atom  abstraction. 


FisneZ.  MP2/3-2lG(*)  structures  of  (a)  CjcTS  for  H-atom  abstraction 
by  AiO  from  CtU,  ZAICH  =  106.6®;  (b)  C  CH3AIOH  adduct,  MlOH 
=  128  8®  ZOAIC  =  116.2®,  ZAICH  =  112.2®,  109.4®(2);  (c)  C,  CH3- 
A1(0)H  adduct,  ZOAIH  =  1 13.7®,  ZHAIC  =  127.2®,  zCAlH  =  1 12.4®, 

1 10.3®  (2);  (d)  CCH3OAIH  adduct,  ZHAIO  =  1 14.7® ,  ZAICX:  *  1 38.2®, 

ZOCH  112.2®,  109.5®(2).  Distances  are  in  10"*®  m. 

Berveand  Pettersson  carried  out  a  partial  CASSCF  optimization 
of  an  assumed  C^v  transition  state  (TS)  but  did  not  obtain  enough 
information  to  predict  rate  coefficients.  In  order  to  estinute 
these,  we  have  analyzed  abstraction  by  conventional  transition- 
state  theory  (CTST),  which  requires  vibrational  and  rotational 
data  to  estimate  the  partition  function  of  the  TS.  Accordingly, 
we  reoptimized  the  complete  TS  structure,  at  the  MP2/3- 
21G{^)  level, 2*  and  obtained  the  vibrational  frequencies.  The 
geometry  is  shown  in  Figure  2a  and  confirms  their*  work;  the 
frequencies  are  1574/,  175(2),  441(2),  570, 1030, 1357(2),  1394, 
1561(2),  3109  and  3226(2)  cm-*.  For  an  energy  barrier  at  0  K 
(i.c.,  including  zero-point  energy,  ZPE)  of  £0,  the  CTST 
abstraction  rate  coefficient  is 

*(59(3-1380  K)  =  1.43  X  lO-^Cr/K)'  '"  x 

cxp(-940  K/T)exp{-Eo/RT)  cm’  molecule"'  s"'  (4) 

Preliminary  estimates  showed  that  quantum  mechanical  H-atom 
tunneling  has  little  effect  when  *  >  1  x  lO""*  cm^  molecule*'  s*'. 
Agreement  between  eqs  2  and  4  at  825  K,  the  midpoint  of  the 
range  of  1^'  investigated,  is  obuined  when  £o(4)  =  1 2.4  kJ  mol*' . 
This  yields  *(  1 380  K)  within  5%  of  the  experimental  value,  while 
*(590  K)  is  too  small  by  a  factor  of  2.3.  Though  accord  at  low 
temperatures  might  be  improved  by  adjustment  of  the  TS 
frequencies,  the  fitted  £o(4)  is  much  lower  than  the  £o(l)  =  67 
kJ  mol*'  from  the  earlier  ab  initio  estimate.'  Use  of  that  value 
together  with  the  observed  total  AIO  removal  rate  coefficient,  eq 
2,  implies  a  branching  ratio  for  direct  abstraction  of  about  9  X 
1 0-’  at  1 380  K  and  about  6  X  10"*  at  590  K.  In  order  to  estimate 
the  uncertainty  in  the  £o(l)  value  that  Borve  and  Pettersson 
derived  at  the  CCI  -I-  Q  level  of  theory,  we  compare  their  estimates 
of  AW*  for  reaction  I  and  the  Li  and  Mg  analogues'  with  the 
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experimental  values.-^  The  mean  error  is  about  ±25  IcJ  moL', 
which  is  a  minimum  estimate  of  the  uncertainly  in  £o(U- 
also  allow  for  a  possible  factor  of  3  uncertainty  in  the  TS  partition 
functions.  These  uncertainties  imply  that  branching  ratios  for 
direct  abstraction,  i.e.,  without  formation  of  an  intermediate  bound 
complex,  as  large  as  0.23  at  1380  K  and  0.003  at  590  K  cannot 
be  excluded,  but  another  reaction  channel,  which  is  totally 
dominant  at  the  lower  temperatures,  is  definitely  indicated. 

The  alternative  abstraction  channel  leading  to  HAIO  +  CH3 
would  be  192  ±  93  kJ  moL‘  endothermic^-^  and  can  on  that  basis 
be  excluded  from  further  consideration.  The  likelier  prospects 
for  further  channels  involve  additions  of  AlO  to  CH4-  Particularly , 
insertion  of  AlO  (a  doublet  species)  into  a  C--H  bond  offers  an 
intriguing  prospect,  by  analogy  to  the  suggested  insertion  of  A1 
and  other  group  1 3  doublet  ground-state  atoms  into  such 
To  test  this,  we  have  employed  theoretical  methods  to  characterize 
possible  insertion  products  and  estimated  the  likely  pressure 
dependence  of  these  to  sec  whether  they  are  consistent  with  our 
observations.  Ab  initio  geometries  of  three  potential  insertion 
species  obtained  by  means  of  MP2/3-21G(*)  theory  are  shown 
in  Figure  2b-d.  Structure  d  was  previously  proposed  as  the 
product  of  insertion  of  A1  into  CH30H,*^  similar  to  «nch  insertion 
into  To  confirm  that  these  structures  are  true  energy 

minima,  we  calculated  the  vibrational  fre<^ucncies,  which  arc  all 
real.'^  These  frequencies  arc  also  used  to  derive  the  zero-point 
energy  (ZPE),  to  obtain  the  temperature  dependence  of  the 
enthalpies  in  the  kinetic  calculations  outlined  below.  The 
enthalpies  of  CH3AIOH,  CH3A1(0)H,  and  CH3OAIH,  relative 
to  AlO  +  CH4  at  0  K,  arc  predicted  at  the  MP4SDQ/6- 
3lG*+ZPE//MP2/3-21G(*)  level  to  be  ^166,  -46,  and  -87  kJ 
moL ' ,  respectively.’^  Equilibrium  constant  calculations  with  these 
enthalpies  indicate  that  CH3AIOH  is  stable  under  our  experi¬ 
mental  conditions,  that  the  second  adduct  is  too  wealdy  bound 
to  be  a  significant  final  sink  for  AlO,  and  that  formation  of  the 
third  adduct  cannot  be  ruled  out  at  the  lowest  temperatures 
studied.  There  is  a  possibility  that  these  adducts  decompose  to 
new  products.  For  example,  CH3AIOH  might  dissociate  to  AlOH 
+  CH3,  perhaps  with  a  lower  overall  energy  barrier  than  proposed 
for  direct  abstraction. 

Next,  the  likely  pressure  dependence  of  AlO  +  CH4 
C  H  3  AlOH  at  8  25  K  was  investigated  by  means  of  QRRK  theory  P 
The  addition  rate  coefficient  is  predicted  to  be  within  25%  of  the 
high-pressure  limit  at  1 5  mbar  of  Ar  and  would  show  only  a  small 
pressure  variation  over  the  experimental  range  with  assumed 
Arrhenius  parameters  for  insertion  of  A  —  10"‘^  cm^  molecule^’ 
S'  ’  and  £a  =  0.  This  is  consistent  with  the  lack  of  an  observed 
pressure  dependence  of  the  experimental  rate  coefficients.  If 
there  were  a  significant  barrier  to  insertion  (£,  >  0),  then 
dissociation  of  the  excited  adduct  back  to  AlO  +  CH4  would  be 
less  competitive  with  collisional  stabilization  and  the  predicted 
rate  coefficients  would  be  even  closer  to  the  high-pressure  limit. 
The  observed  pressure  independence  is  also  consistent  with  a 
mechanism  where  the  adduct  fragments  to  new  producu  before 
stabilization.  Further  QRRK  calculations  suggest  that  formation 
of  any  adduct,  which  is  sufficiently  thermodynamically  stable  to 
be  significant  in  our  experiments,  will  have  a  rate  coefficient  that 
is  in  the  falloff  region  or  close  to  the  high-pressure  limit. 

Cooclusions 

These  experiments  arc  in  accord  with  the  suggestion  by  Borve 
and  Pettersson*  that  gas-phase  reactions  between  AlO  and  CH4 
occur  and  may  provide  insight  in  the  heterogeneous  catalytic 
conversion  of  methane  to  C2  hydrocarbons.  However,  use  of 
their  calculated  barrier  suggests  that  the  branching  ratio  for  direct 
abstraction  is  small  and  therefore  that  other  reaction  channels 
are  important.  Our  calculations  demonstrate  that  insertion 
leading  to  CH3AIOH  is  a  reasonable  candidate  for  the  iniual 
step,  although  other  channels  arc  not  excluded.  The  insertion 
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products  (Figure  2b-d)  might  model  the  products  of  dissociative 
adsorption  of  CH4  onto  an  alumina  catalyst,  where  their  energies 
arc  likely  to  differ  from  gas-phase  values.  More  detailed 
experimental  and  theoretical  investigations  of  the  present  reaction 
system,  as  well  as  of  other  metal  oxide-hydrocarbon  rcaaions, 
arc  needed  to  obtain  an  understanding  of  this  new  class  of  gas- 
phase  reactions  and  their  potential  importance  for  methane 
conversion. 
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The  high-temperature  fast-flow  reactor  technique  has  been  used  to  make  kinetic  measurements.  A  weighted 
fit  to  the  AlCl  +  N2O  data  gives  Jt(700-990  K)  =  5.6  X  10-»>  exp(-7380  K/T)  cm^  molecule-'  s-'.  A  weighted 
fit  to  the  AlCl  -h  CO2  measurements  leads  to  the  expression  *(900-1790  K)  =  4.4  X  10-^3  (F/K)  •  exp(-3^ 
K/D  cm^  molecule-'  s-'.  2<r  accuracy  limits  are  about  ±25%.  An  upper  limit  *(800-1100  K)  <  5  X  10- 
cm^  molecule-'  s-'  has  been  determined  for  the  AlCl  +  SO2  reaction.  An  alternate  form  of  classical  transition- 
state  theory  is  developed  to  allow  predictions  on  the  preexponential  part  of  rate  coefficient  expressions  for 
metallic  spwies.  This  model-based  transition-state-theory  (MTST)  method  uses  a  valence-force  molwular 
model  to  estimate  rotational  constants  and  vibrational  frequencies  of  the  transition  sute  and  is  applicable  to 
reactions  with  early  barriers,  typical  of  many  exothermic  charge-transfer  reactions.  The  grometncal  paramet^ 
and  force  constants  that  describe  the  molecular  model  are  derived  from  properues  of  the  reactants.  For  the 
NiO  reaction  good  agreement  between  MTST  and  experiment  is  obtained,  based  on  the  assumption  of  an  O 
atom  abstraction  reaction  leading  to  OAlCl.  No  such  agreement  is  found  for  the  CO2  indiMtes 

adduct  formation  as  the  main  AlCl  consumption  channel.  For  the  previously  measur^  AlCl  +  O2 
MTST  calculations  suggest  that  abstraction  can  be  of  some  significance  above  about  1 500  K;  however,  adduct 
formation  appears  to  dominate  over  most  of  the  490-1750  K  range. 


Introduction 

As  part  of  an  ongoing  study  of  the  temperature  dependence 
of  the  rate  coefficients  of  reaaions  of  small  Al-  and  B-containing 
species,  we  report  here  experiments  on  the  reactions 


AlCl  +  NjO  -*  products 

(1) 

AlCl  +  CO2  -*  products 

(2) 

AlCl  +  SO2  —  products 

(3) 

We  previously  studied  the  AlCl  reactions  with  02,^  C02,^  Cl2i^ 
and  HCl*  the  results  of  which  have  been  summarized.^  The 
measurements  of  reaction  2  are  repeated  here,  as  the  high- 
temperature  fast-flow  reactor  (HTFFR)  technique  used  has  been 
much  improved^’^  since  the  previous^  study,  where  the  data  showed 
wide  scatter.  The  O2  and  CO2  reactions  were  originally  thought 
to  lead  direaly  to  OAlCl.  However,  a  recent  ab  initio  study  of 
the  thermochemistry’  of  OAlCl  suggests  these  two  channels  to 
be  too  endothermic*  for  reaction  along  this  path  with  the  observed 
rate  coefficients. 

To  investigate  the  mechanism  of  all  these  AlCl  reactions  with 
oxygen-containing  oxidants,  we  present  here  a  semiempirical 
transition-state  theory  based  method  for  predicting  the  temper¬ 
ature  dependence  of  rate  coefficients.  In  this  model-based 
transition-state  theory  (MTST)  the  rotational  constants  and 
vibrational  frequencies  of  the  transition  state  arc  derived  from 
molecular  models  based  on  the  assumption  of  valence  forces. 
This  approach,  the  limitations  of  which  are  discussed,  should  be 
widely  applicable  to  reaaions  with  early  activation  barriers  since 
the  required  input  data,  namely,  geometrical  parameters  and 
force  constants,  of  the  transition  state  can  be  estimated  from 
reactant  properties  in  accord  with  Hammond’s  postulate.^ 
Reactions  occurring  by  a  crossing  of  neutral  and  ionic  potential 

^  Present  address:  Department  of  Chemical  Engineering,  Massachusetts 
Institute  of  Technology,  Cambridge,  MA  02139. 
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energy  surfaces,  typical  of  oxidation  reactions  of  metal  species, 
are  included.  MTST  is  thus  useful  for  such  reaaions  when  the 
input  data  for  the  widely  used  semiempirical  TST  method 
developed  by  Benson,^®  principally  for  H/C/O/N  species  reac¬ 
tions,  arc  not  available. 

Experimental  Section 

Detailed  descriptions  of  the  HTFFR  used  in  this  work  have 
been  given  elsewhere.^*^  The  reactor  consists  of  a  vertical  muUitc 
reaction  tube  (2.2-cm  i.d.)  surrounded  by  SiC  heating  elements. 
This  assembly  is  sealed  inside  an  insulated  water-cooled  vacuum 
housing.  Two  methods  were  used  to  produce  AlCl.  In  the  first, 
previously  used  method,  a  trace  amount  of  CI2  or  HCl  was  added 
to  the  Ar  bath  gas  flowing  over  the  Al-wetted  resistively-heated 
tungsten  vaporization  coil,  located  7—17  cm  upstream  from  the 
reaction  zone.^’^  In  the  second  method,  similar  to  that  used  in 
the  BCl  studies,^»^^  AICI3  was  vaporized  from  a  quartz  crucible 
(4-mm  i.d.)  into  a  stream  of  Ar.  This  stream  then  flowed  through 
a  2450’MHz,  nominally  1 00- W,  microwave  discharge  in  a  quartz 
tube  ( 1 .2-cm  i.d.)  and  then  into  the  reaction  tube.  Both  methods 
were  used  in  the  AlCl  +  CO2  study  in  the  1200-1690  K 
temperature  range  and  found  to  yield  identical  results.  The  second 
method  would  in  general  be  especially  useful  for  HTFFR  studies 
near  room  temperature;  such  temperatures  are  not  accessible  by 
the  first  method  because  the  heat  generated  by  the  vaporization 
coil  heats  the  reaaion  tube  to  temperatures  greater  than  450  K. 
Downstream  from  the  source,  a  mixture  of  Ar  and  the  oxidant 
gas  (N2O,  CO2,  or  SO2)  was  introduced  into  the  reaction  tube 
through  a  movable  inlet.  The  molar  flow  rate  of  gas  through  the 
inlet  was  kept  constant  during  a  given  rate  coefficient  determi¬ 
nation,  which  consisted  of  five  or  six  measurements  at  varying 
oxidant  concentrations.  A  small  amount  of  O2  was  added  to  the 
inlet  gas  to  rapidly^  remove  any  remaining  Al  atoms  entering  the 
reaction  zone.  Such  atoms  otherwise  may  continue  to  produce 
AlCl  and  hence  interfere  with  the  measurements. 

Relative  AlCl  concentrations  were  monitored  by  laser-induced 
fluorescence  using  a  pulsed  Lambda  Physik  EMG  101  excimer/ 
FI  2002  dye  laser  combined  with  a  KDP  doubling  crystal.  The 
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TABLE  I:  Summary  of  Rate  Coefficient  Measurements  for 
AlCl  +  N2O-  ' 


reaction 


zone 

length 

(cm) 

P 

(mbar) 

(Ml 

(10‘’ 

cm~^) 

[N20]m« 

(101^ 

cm"^) 

F 

(m  H) 

T 

(K) 

k±  Ck  (cm^ 
molecule-*  s"*) 

20 

50.9 

4.0 

16.4 

40^ 

20 

919 

2.01  ±0.13  (-14)^ 

10 

50.9 

4.0 

16-5 

73 

20 

919 

2.16  ±0.17  (-14) 

20 

50.9 

4.0 

13.7 

36 

21 

921 

2.01  ±0.12  (-14) 

10 

50.9 

4.0 

13.6 

27 

21 

924 

2.18  ±0.19  (-14) 

10 

101.0 

7.7 

14.1 

34 

14 

948 

2.29  ±0.22  (-14) 

10 

101.2 

7.6 

13.9 

37 

14 

964 

2.08  ±0.22  (-14) 

20 

29.6 

2.2 

14.3 

50 

39 

984 

2.75  ±0.21  (-14) 

20 

29.7 

2.2 

14.2 

62 

39 

991 

2,82  ±0.25  (-14) 

20 

54.1 

5.3 

23.9 

23 

16 

735 

2.40  ±0.33  (-15) 

20 

54.1 

5.4 

24.2 

22 

16 

726 

2.19  ±0.32  (-15) 

20 

78.6 

8.1 

52.3 

30 

7 

704 

1.44  ±0-14  (-15) 

20 

78,6 

8.1 

52.3 

17 

7 

704 

1.80  ±0.13  (-15) 

20 

66.7 

6.9 

40.0 

49 

14 

704 

1.31  ±0.22  (-15) 

20 

66.7 

6.8 

39.7 

48 

14 

709 

1.11  ±0.22  (-15) 

20 

41.6 

4.2 

22.2 

35 

13 

726 

2-23  ±0.22  (-15) 

20 

41.6 

4.1 

22.1 

33 

13 

731 

2.45  ±0.21  (-15) 

20 

55.2 

5.4 

28.9 

34 

10 

740 

2.80  ±0.21  (-15) 

20 

55.2 

5.3 

28.6 

41 

10 

748 

3.13  ±0.22  (-15) 

20 

78.5 

7.2 

35.6 

31 

12 

795 

4.17  ±0,35  (-15) 

20 

62.9 

5.5 

32.7 

50 

17 

830 

7.68  ±0.72  (-15) 

20 

62.9 

5.4 

32.2 

54 

17 

842 

.9.51  ±0.64  (-15) 

10 

62.9 

5.5 

32.3 

59 

17 

836 

9.35  ±0.67  (-15) 

10 

49.9 

4.2 

15.0 

23 

19 

856 

1.04  ±0.08  (-14) 

20 

49.9 

4.2 

15.1 

26 

19 

853 

9.80  ±0.70  (-15) 

20 

42.5 

3.6 

13.5 

31 

21 

846 

9.20  ±0.79  (-15) 

20 

42.8 

3.6 

13.5 

24 

21 

852 

1.08  ±0.08  (-14) 

20 

44.3 

3.6 

17.3 

42 

32 

879 

1.30  ±0.16  (-14) 

20 

44.3 

3.6 

17.3 

34 

32 

882 

1.48  ±0.10  (-14) 

20 

35.2 

2.9 

13.9 

17 

40 

873 

1.50  ±0.15  (-14) 

20 

50.7 

4.2 

20.0 

9 

27 

871 

8.86  ±0.82  (-15) 

20 

68.4 

5.7 

17.3 

52 

16 

875 

1.07  ±0.07  (-14) 

®  The  measurements  arc  reported  in  the  sequence  in  which  they  were 
obtained.  *  In  arbitrary  units.  ^  Should  be  read  as  (2.01  ±  0.13)  x  10"^^ 
cm^  molecule'^  s“^ 


AlCUA^n-X^S*^)  transitions^  was  pumped  on  the  (0,0)  band  at 
261.4  nm,  and  the  fluorescence  intensity  was  observed  through 
a  262-nm  (26-nm  fwhm)  filter  by  an  EMI  9813QA  photomul¬ 
tiplier  tube  coupled  to  a  Data  Precision  Analogic  6000/620  100- 
MHz  transient  digitizer.  Rate  coefficient  measurements  were 
made  in  the  stationary  inlet  mode^^  with  observed  reaction  zone 
lengths  of  10  or  20  cm.  The  oxidant  concentrations  were 
maintained  in  great  excess  over  the  AlCl  concentrations.  Rate 
coefficients  and  their  uncertainties  were  obtained  from  the  slopes 
of  linear  plots  of  In  [AlCl] relative  vs  oxidant  concentration  by  a 
weighted  linear  regression. 

The  gases  used  were  99.998%  Ar  from  the  liquid,  “Precision 
Aquarator”  99.99%  CO2,  CI2  (1%  in  Ar),  and  “Atomic  Absorp¬ 
tion”  99.0%  N2O  from  Linde,  “Ultra  High  Purity”  99.99%  NjO 
and  “Anhydrous”  99.98%  SO2  from  Matheson,  HCl  (5%  in  He) 
from  Sp)ectra,  and  O2  (5.0%  in  Ar)  from  Scott. 

Results 

The  AIQ  +  N2O  Reaction.  Reaction  1  was  studied  from  700 
to  990  K.  No  measurements  at  higher  temperatures  are  included 
as  thermal  dissociation  of  N2O  is  sufficiently  large  to  interfere 
with  rate  coefficient  measurements. •  ‘  Below  700  K  the  reaction 
was  too  slow  to  measure.  Table  I  summarizes  the  rate  coefficients 
and  the  temperatures  at  which  they  were  obtained  along  with  the 
experimental  parameters:  reaction  zone  length,  total  pressure  P, 
total  concentration  [M],  maximum  oxidant  concentration 
[N20]ni«»  fluorescence  intensity  F  (a  measure  of  [AlCl]  at  the 
start  of  the  reaction  zone),  and  average  velocity  0.  The 
independence  of  the  individual  rate  coefficients  A:,  of  these 
parameters  was  checked  by  examining  plots  of  [k{T)-  ki]/k(T), 
where  k(  T)  is  obtained  from  the  fitted  rate  coefficient  expression 
given  below,  vs  the  above  experimental  parameters.  It  can  be 
seen  from  Figure  1  that  over  the  limited  temperature  range  studied 


T{K) 

1300  1000  700 


Figure  1.  Arrhenius  plots  of  cxpcrimcnul  rate  coefficients  for  the  AlCl 
+  N2O  (open  circles)  and  AlC!  +  CO2  (open  squares)  reactions.  The 
experimental  rate  coefficient  expressions  arc  shown  as  full  lines, 

the  ki  measurements  show  no  deviation  from  Arrhenius  behavior. 
A  weighted  nonlinear  fit^^  of  the  data  to  the  expression  A(T)  = 
A  expi-E/RT)  gives 

yt, (700-990  K)  =  5.6  x  10"' ‘  X 

exp(-7380  K/T)  cm^  molecule”’ s‘‘  (4) 
with  covariance  matrix  elements^® 

{ajAf  =  Q.\U,  ^^/  =  79  900K^  a^lA  =  9^.€Yi 

(5) 

Combining  these  matrix  elements  by  the  propagation  of  errors 
technique^-^^  yields  ±2(7*  precision  limits  of  15%  at  700  K,  6% 
at  850  K,  and  1 2%  at  990  K.  Allowing  10%  uncertainty  associated 
with  the  flow  profile  factor  and  20%  for  further  systematic  errors 
yields  ±2cr*  accuracy  limits  of  27%  at  700  K,  23%  at  850,  and 
25%  at  990  K. 

The  AlCl  +  CO2  Reaction.  Reaction  2  was  studied  over  the 
900-1790  K  temperature  range.  The  rate  coefficient  data 
summarized  in  Table  II  are  again  independent  of  reaction  zone 
length,  P,  [M],  (C02]ni«,  F,  and  v.  The  Arrhenius  plot  of  the 
rate  coefficients  shown  in  Figure  1  exhibits  slight  upward 
curvature.  The  rate  coefficient  data  were  fitted  to  the  three- 
parameter  expression  k(T)  -  AT^  exp(-£/P71.  The  goodness 
of  the  fit  is  found  to  be  relatively  insensitive  to  the  value  of 
the  parameter  n  over  the  approximate  range  2.0-4.0.  We  therefore 
fix  n  at  3.0.  A  two- parameter  nonlinear  regression  analysis  of 
the  data  gives 

^2(900-1790  K)  =  4.4  X  10”^^  (r/K)^  “  X 

exp(-3900  K/  T)  cm^  molecule”’  s”’  (6) 
with  covariance  matrix  elements 

((rV/4)^  =  0.0182,  (r£^  =  38  000K^  =  26.0  K 

(7) 

Combining  the  covariance  matrix  elements  by  the  propagation 
of  errors  technique  yields  ±2(t*  precision  limits  of  17%  at  900  K, 
5%  at  1350  K,  and  6%  at  1790  K.  Allowing  10%  uncertainty 
associated  with  the  flow  profile  factor  and  a  conservative  20%  for 
further  systematic  errors  yields  ±2(7*  confidence  limits  of  28% 
at  900  K,  23%  at  1350,  and  23%  at  1790  K.  The  present  rate 
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TABLE  11:  Summary  of  Rate  Coefficient  Measurements  for 
AlO  +  CO2'  ^  _ 


reaaion 


zone 

length 

(cm) 

P 

(mbar) 

cm-’) 

[C02]„„ 

(10^5 

cm-3)  F 

V 

(m 

s-‘) 

T 

(K) 

ck  (cm^ 
molecule"*  s"*) 

20 

43.3 

2.5 

18 

18* 

25 

1273 

4.20 

±0.55(-15)‘^ 

20 

78.4 

4.5 

33 

71 

13 

1254 

3.73 

±0.50  (-15)-' 

20 

78.6 

4.5 

32 

64 

14 

1273 

3.21 

±0.22  (-15)^ 

10 

78.6 

4.4 

32 

62 

14 

1281 

3.86 

±0.53  (-15)-' 

10 

78.6 

4.4 

32 

108 

14 

1284 

4.15 

±0.28  (-15)-' 

10 

63.5 

3.2 

17 

33 

14 

1429 

8.71 

±0.98  (-15)“' 

10 

63.7 

3.2 

21 

47 

14 

1440 

9.99 

±0.91  (-15)“' 

20 

62.8 

3.1 

12 

50 

14 

1475 

9.69 

±0.98  (-15)“' 

20 

62.8 

3.0 

12 

54 

14 

1493 

1.28 

±0.13  (-UV* 

20 

43.5 

2.8 

16 

13 

19 

1132 

2.84 

±0.30  (-15)-' 

20 

43.5 

2.7 

15 

18 

19 

1149 

2.38 

±0.98  (-15)^ 

20 

58.1 

4.3 

23 

10 

13 

990 

7.83 

±  2.20  (-16H 

20 

57.7 

4.1 

23 

23 

13 

1010 

1.06 

±0.18  (-15)^ 

20 

88.1 

6.2 

27 

27 

11 

1031 

1.15 

±0.17  (-15)*' 

20 

46.4 

3.1 

20 

59 

22 

1091 

1.27 

±0.34  (-15)'* 

20 

46.1 

3.1 

20 

48 

22 

1094 

1.34 

±0.45  (-15)*' 

20 

53.5 

2.4 

14 

21 

31 

1584 

1.33 

±  0.19  (-MH 

10 

53.3 

2.4 

14 

24 

31 

1602 

1.74 

±  0.21  (-14)^ 

10 

53.2 

2.4 

14 

21 

31 

1609 

2.04 

±  0.18  {-I4y 

10 

92.0 

4.0 

20 

21 

15 

1679 

1.0/ 

±  0.24  (-14)^ 

10 

92.0 

4.0 

20 

17 

15 

1686 

2.05 

±  0.21  (-14)^ 

20 

53.7 

4.2 

21 

55 

11 

930 

7.20 

±  3.38  (-16)<' 

20 

53.7 

4.1 

21 

43 

11 

942 

9.32 

±3.11  {-ley 

20 

65.3 

5.1 

27 

74 

16 

920 

5.70 

±2.62  (-16)“' 

20 

65.7 

5.1 

27 

54 

16 

925 

6.69 

±  1.47  (-16)'' 

20 

78.6 

6.3 

36 

39 

16 

900 

5.97 

±  1.34  (-16H 

20 

43.1 

2.3 

17 

19 

34 

1368 

5.49 

±  1.38  (-15)^ 

20 

43.1 

2.3 

17 

40 

35 

1373 

5.24 

±0.71  (-15)-' 

10 

43.1 

2.3 

17 

68 

35 

1374 

6.10 

±2.76  (-15)-' 

10 

43.1 

2.3 

17 

66 

35 

1386 

6.77 

±0.83  (-15)-' 

10 

85.8 

4.5 

33 

40 

18 

1397 

6.66 

±0.68  (-15)-' 

10 

85.4 

4.4 

33 

19 

18 

1399 

6.00 

±0.95  (-15)^ 

20 

83.6 

4.3  ■ 

17 

30 

18 

1401 

5.74 

±  1.00  (-15)-' 

20 

99.6 

6.0 

30 

113 

18 

1210 

3.69 

±0.41  (-15)' 

20 

113.3 

6.8 

30 

87 

18 

1204 

2.72 

±0.28  (-15)' 

20 

73.3 

4.2 

25 

59 

17 

1250 

3.75 

±0.32  (-15)' 

20 

60.4 

2.8 

12 

36 

45 

1540 

1.81 

±0.15  (-14)' 

20 

60.4 

2.8 

12 

39 

46 

1556 

1.80 

±0.23  (-14)' 

20 

91.4 

4.2 

16 

39 

25 

1590 

1.83 

±0.10  (-14)' 

20 

91.6 

4.2 

16 

78 

25 

1590 

1.90 

±0.11  (-14)' 

10 

92.8 

4.2 

21 

53 

25 

1599 

1.71 

±  0.12  (-14)' 

10 

93.0 

4.2 

21 

54 

25 

1599 

1.74 

±0.12  (-14)' 

10 

61.9 

2.7 

18 

62 

22 

1632 

2.14 

±0.15  (-14)' 

20 

60.5 

2.7 

12 

65 

22 

1638 

2.32 

±0.14  (-14)' 

20 

52.9 

2.6 

17 

87 

33 

1448 

6.59 

±0.54  (-15)' 

20 

52.8 

2.6 

16 

59 

33 

1449 

6.98 

±0.64  (-15)' 

10 

52.8 

2.6 

16 

60 

33 

1454 

9.33 

±  1.11  (-15)' 

10 

52.8 

2.6 

16 

60 

33 

1454 

9.08 

±  1.06  (-15)' 

10 

75.0 

3.1 

12 

70 

37 

1765 

2.75 

±0.28  (-14)' 

10 

75.2 

3.1 

12 

69 

36 

1765 

2.81 

±0.20  (-14)' 

20 

74.8 

3.0 

7 

87 

37 

1777 

2.26 

±0.14  (-14)' 

20 

74.8 

3.0 

12 

88 

37 

1777 

2.57 

±0.24  (-14)' 

20 

41.9 

1.7 

8 

91 

68 

1787 

2.61 

±  0.21  (-14)' 

20 

42.1 

1.7 

8 

125 

67 

1787 

2.55 

±0.16  (-14)' 

20 

42.1 

1.7 

8 

39 

67 

1791 

2.98 

±0.26  (-14)' 

^  The  measurements  are  reported  in  the  sequence  in  which  they  were 
obtained.  *  In  arbitrary  units.  ^  Should  be  read  as  (4.20  ±  0.55)  x 
cm^  molecule*^  s"‘ .  ^  AlCl  produced  by  the  microwave-discharge  method. 
'  AlCI  produced  by  the  reacting  CI2  or  HCl  with  an  Al-wetted  vaporization 
coil. 


coefficient  recommendation  is  in  good  agreement  with  that 
determined  in  the  previous  investigation;^  however,  the  scatter  in 
the  measurements  is  now  greatly  reduced. 

The  AIQ  +  SO2  Reaction.  The  rate  coefficients  of  reaction 
3  could  not  be  accurately  determined  because  of  absorption  and 
emission  of  laser  radiation  by  S02,^^  which  greatly  lowered  the 
sensitivity  of  AlCl  detection.  Corrections  were  made  to  the 
iniegratwi  fluorescence  signal  to  remove  the  SO2  contribution. 
Over  the  800-1 100  K  temperature  range  no  significant  change 
in  the  AlCl  signal  could  be  observed  as  the  SO2  concentration 
was  varied,  which  indicates  an  upper  limit  to  the  rate  coefficient 
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jfe3(800-1100K)<5xiO'“‘  cm^  molecule-*  s-‘  (8) 

We  were  unable  to  investigate  reaction  3  at  higher  temperatures 
because  the  emitted  radiation  from  species  other  than  AlCl  was 
too  intense.  Thermochemical®  equilibrium  calculations  indicate 
that  above  1 100  K  SO,  which  also  absorbs  and  emits  radiation 
near  the  AlCl  transition,  5^  present  in  significant  quantities. 

Model-Based  Transition-State  Theory  (MTST) 

According  to  conventional  transition-state  theory  the  rate 
coefficient  for  a  reaction  occurring  on  a  single  potential  energy 
surface  is  given  by^’ 

(9) 

in  which  qx  is  the  partition  function  for  the  transition  state,  q*A 
and  qb  are  partition  functions  for  the  reactants,  and  is  the 
difference  in  zero-point  levels  of  the  transition  state  and  the 
reacunts,  and  are  readily  evaluated  from  experimental 
values  of  the  rotational  constants  and  vibrational  frequencies  of 
the  reactants.  The  data  needed  to  evaluate  qx  are  seldom  available 
and  must  therefore  be  estimated.  The  MTST  approach  for  this 
evaluation  is  suitable  for  transition  states  which  arise  from  an 
early  crossing  of  neutral  and  ionic  potential  surfaces,  such  as  is 
the  case  for  many  metal  species  reactions.  For  such  a  transition 
state  reasonably  good  estimates  of  the  geometrical  parameters 
and  vibrational  frequencies  can  be  made  using  Hammond's 
postulate.^  It  can  then  be  assumed  that  the  intemuclear  distances, 
valence  angles,  and  force  constants  of  the  transition  states  arc 
closely  related  to  the  properties  of  the  reactants  for  an  exothermic 
reaction  and  properties  of  the  products  for  an  endothermic 
reaction. 

Vibrational  frequencies  for  transition  states  are  obtained  from 
molecular  models.  The  molecular  model  used  here  assumes  that 
harmonic  restoring  forces  oppose  changes  in  bond  lengths,  angles 
between  valence  bonds,  and  dihedral  angles.  The  potential  energy 
V  in  the  region  near  the  transition  state  is  thus  given  by*® 

2K  =  (10) 

1  j  k 

where  Fb,  and  Fr  are  force  constants  for  changes  in  intemuclear 
distances  r,  valence  angles  and  dihedral  angles  a  from  their 
respective  equilibrium  values.  A  normal-coordinate  analysis  of 
the  model,  which  procedure  has  been  described  in  detail  by  Pauling 
and  Wilson,*^  is  carried  out.  The  secular  equation  thus  obtained 
is  solved  numerically  by  reducing  the  determinant  to  the 
tridiagonal  form  using  Householder’s  algorithm.*^  The  char¬ 
acteristic  polynomial  of  the  resulting  matrix  is  then  evaluated  by 
the  QL  algorithm  with  implicit  shifts.*^  Vibrational  frequencies 
Vi  are  obtained  from  the  relation*’ 

;.,=  X“-'/2tc  (11) 

in  which  X,  are  the  eigenvalues  of  the  secular  equation.  Rotational 
constants  B,-  for  the  transition  state  are  obtained  from  the 
eigenvalues  of  the  inertia  tensor,^  that  is,  from  the  principal 
moments  of  inertia  /,  by  the  relation^* 

5,  =  A/(8x^c/,)  (12) 

Several  examples  of  accordingly  estimating  transition-state 
parameters  for  AlCl  reactions  are  given  in  the  Discussion. 

Some  limitations  of  MTST  are  that  force  constants  arising 
from  the  reaction,  which  includes  one  stretching  force  constant 
and  one  or  more  bending  or  torsional  force  constants,  must  be 
estimated  from  the  available  ab  initio  calculations  on  similar 
reactions.  The  most  important  of  these  are  the  torsional  force 
constants,  since  the  torsions  can  be  of  low  frequency  and  hence 
influence  the  TST  rate  expression.  Also,  MTST  is  not  suitable 
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TABLE  III:  Transirion-State  Parameters  for  ClAJONN* 

geometricai  parameter®'^  force  constant‘s 


r(ClAl) 

2.13 

F(CIAI) 

2.1 

r(AlO) 

2.25 

F.(A10) 

-2.0 

r(ON) 

1.18 

Fr(ON) 

5.85 

r(NN) 

1.13 

F.(NN) 

18.7 

^(CIAIO) 

120 

Fs(ClAlO) 

0.2 

^(AlON) 

90 

F^(AION) 

0.6 

^(ONN) 

160 

FflCONN) 

0.67 

r(ClAiON) 

0 

F,(C1A10N) 

0.02 

r(AlONN) 

0 

F,(AIONN) 

0.02 

rotational  constants  {B1B2BJ):  9.51  X  1(H  cm  ^ 
vibrational  frequencies  Vi:  39,  89,  150,  291,  478,  605,  991,  2221  cm“^ 

^  We  assume  trans  arrangements  for  the  atoms  in  the  CIAION  and 
.AlONN  groups.  *  Units  of  r  are  10“^®  m;  units  of  6  and  r  arc  deg.  ‘s  Units 
of  Fr  are  10^  N  m"^;  units  of  F$  and  Fr  arc  10“**  N  m  rad“2. 

for  reactions  with  a  centrally  located  transition  state,  e.g.,  H  + 
HCl  —  H2  +  Cl,  because  the  potential  energy  function  (cf.  eq 
10)  would  have  to  be  modified  to  include  a  cross-term  for  the 
interaction  between  the  H-H  and  H-Cl  bonds.^^  Furthermore, 
the  transition  state  would  resemble  both  the  reac  ints  and  the 
products  in  part;  hence,  there  would  be  considerable  uncertainty 
in  estimating  geometrical  parameters  and  force  constants. 

Discussion 

The  AIG  +  N2O  Reaction.  The  lowest  channel  of  the 
abstraction  reaction 

AlCl  +  N.O  —  O  AlCl  +  N2  (la) 

is  spin-allowed  and  267  kJ  moh^  exothermic.-^  The  parameters 
of  the  transition  state  CIAIONN^  are  derived  from  properties  of 
the  reactants  and  are  summarized  in  Table  III.  The  distances 
r(ClAl),  r(ON),  and  r(NN)  are  from  values  for  the  AlCl  and 
N2O  molecules.  The  internuclear  distance  along  the  reaction 
coordinate  r(AlO)  is  estimated  by 

r(AlO)  =  [KAD  +  r(0-)]0.50‘/^^'’  (13) 

Equation  1 3  is  based  on  a  relation  given  by  Pauling^^  between 
univalent  radii  and  the  bond  distances  of  diatomic  alkali-metal 
halide  molecules.  We  assume  that  the  relation  applies  to  charge- 
transfer  complexes  and  hence  to  transition  states,  which  are  located 
fairly  close  to  the  complexes.  Values  for  the  univalent  radii  KA^) 
and  r(0“)  and  the  average  Bom  exponent  5  are  obtained  as 
described  in  ref  25,  The  constant  0.50  is  obtained  such  that  eq 
1 3  gives  good  agreement  with  ab  initio  computations^^^  of  B-0, 
Al^,  and  H-O  distances  for  transition  states  of  the  respective 
BCl  +  O2,  A1  -h  CO2,  and  H  -H  N2O  reactions.  Lacking 
information  on  the  valence  and  dihedral  angles  of  the  transition 
state,  we  assume  the  values  given  in  Table  III  based  on  optimized 
geometries  from  ab  initio  computations. The  force  constant 
F^ClAl)  is  taken  as  that  in  AlCl,^  while  F^(NN)  and  f’^(ONN) 
arc  obtained  from  data  for  N20.^  The  remaining  force  constants 
in  Table  III  are  estimated  by  a  scheme  that  gives  good  agreement 
with  ab  initio  computations^^^  of  vibrational  frequencies.  For 
example,  F^ON)  is  taken  as  one-half  the  N-O  stretching  force 
constant  in  N2O.  In  addition,  all  torsional  force  constants  are 
taken  as  0.02  x  10"**  N  m  rad"^  and  the  stretching  force  constants 
corresponding  to  reaction  coordinates  are  taken  as  ~2.0  X  10^  N 
m"*.  These  data  arc  used  to  calculate  the  transition-state 
rotational  constants  and  vibrational  frequencies,  which  are 
summarized  in  Table  III. 

To  allow  comparison  between  theory  and  experiment,  the 
preexponential  factor  of  eq  9  is  calculated  for  the  500-2000  K 
temperature  range  and  fitted  to  the  expression  A  7^.  To  calculate 
the  barrier  height,  the  experimental  rate  coefficient  expression 
is  set  equal  to  the  TST  equation  (9)  at  820  K,  the  temperature 
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Figure  X  Arrhenius  plots  of  experimental  rate  coefficients  for  the  AlCl 
4“  N2O  (open  circles)  and  AlCl  +  CO:  (open  squares)  reactions.  Rate 
coefficient  expressions  for  the  O  atom  abstraction  channels  for  these 
reactions  derived  using  the  MTST  method  are  shown  as  solid  lines. 

TABLE  rV:  Transition-State  Parameters  for  CIAIOCO* 

geometrical  parameter^  force  constant‘s 


r(ClAl) 

2.08 

F,(C1A1) 

3.2 

r(AlO) 

1.63 

F.(AlO) 

4.2 

riOC) 

1.81 

F.(OC) 

-2.0 

riCO) 

1.13 

FACO) 

19.0 

0(CIA1O) 

160 

F,(CIA10) 

0.3 

<3(AlOC) 

90 

F,(A10C) 

0.6 

6{OCO) 

120 

FeiOCO) 

0.2 

t(CIAIOC) 

0 

FriClMOC) 

0.02 

t(AIOCO) 

0 

Fr(  AlOCO) 

0.02 

rotational  constants  (B1F2B3):  9.46  X  10-*  cm"^ 
vibrational  frequencies  »/,•:  75,  91,  141,  216,  231,  523,  908,  2164  cm"* 

We  assume  trans  arrangements  for  the  atoms  in  the  ClAiOC  and 
AlOCO  groups.  The  parameters  of  the  transition  state  are  based  on  the 
products  of  the  endothermic  reaction  AlCl  +  CO2  CO  +  OAiCl. 
*  Units  of  r  arc  10"*°  m:  units  of  6  and  r  are  deg.  Units  of  Fr  arc  10^ 
N  m"*;  units  of  F0  and  Fr  arc  10"**  N  m  rad"^. 

corresponding  to  the  middle  of  the  T"*  range  of  the  experiments. 
The  value  A£^  =  45.9  kJ  mol"^  is  thus  obtained,  which  yields 

JtMTST  =  5  X  10-'*  X 

exp(-5520  K/T)  cm*  molecule"*  s"'  (14) 

A  comparison  between  eq  14,  the  MTST  rate  coefficient  for 
reaction  la,  and  experiment  is  shown  in  Figure  2.  Over  the 
temperature  range  of  the  experimental  data  the  agreement  can 
be  seen  to  be  excellent,  which  suggests  that  the  dominant  channel 
for  reaction  1  is  abstraction,  i.e.,  channel  la. 

The  AlCl  +  CO2  Reaction.  The  abstraction  reaction 

AlCl  +  CO2  OAlCl  +  CO  (2a) 

is  98  kJ  mol"*  endothermic.^^  The  parameters  of  the  transition 
state  CIAIOCO^  which  are  summarized  in  Table  IV,  have  been 
derived  based  on  properties  of  the  products  OAlCl  +  CO.  The 
estimation  scheme  is  very  similar  to  that  for  the  AlCl  +  N2O 
transition  state  so  only  the  differences  will  be  discussed  here.  The 
r(CO),  KClAl),  and  r(AlO)  distances  are  taken  as  those  in  CO 
and  OAlCl.  The  internuclear  distance  along  the  reaction 
coordinate  is  estimated  by  r(OC)  =  [r(0")  +  r(C‘*")]0.50*/(*"*>. 


7226  The  Journal  of  Physical  Chemistry,  Vol.  97,  No.  28,  1993 

The  force  constant  F,(CO)  is  taken  as  that  in  the  CO  molecule, 
while  FfClAl)  and  F^CCIAIO)  are  taken  as  equal  to  force  constants 
derived  from  ab  initio  calculations'  of  the  vibrational  frequencies 
of  OAlCl.  Table  IV  summarizes  the  resulting  rotational  constants 
and  vibrational  frequencies. 

We  calculate  the  preexponential  factor  of  eq  9  in  the  500- 
2000  K  temperature  range  and  fit  the  results  to  the  expression 
A  expression.  A  barrier  height  of  1 1 8  kJ  moh^  for  the  forward 
reaction  is  estimated  based  on  the  endothermicity  of  reaction  2a 
plus  a  barrier  height  for  the  reverse  reaction  of  20  kJ  mot\ 
which  is  probably  too  small  in  view  of  the  fact  that  CO  +  OAlCl 
is  isoelect ronic  with  AlCl  +  N2O,  for  which  we  obtained  a  barrier 
height  of  45.9  kJ  mol"^ .  Combining  the  calculated  A,  n,  and 
values  leads  to  the  expression 

Atmtst  =  1-27  X  10-'*  {TIKf**  X 

exp(-14190  K/f)  cm*  molecule"' s"'  (15) 

Equation  15  is  plotted  in  Figure  2  to  allow  for  comparison  with 
experiment.  It  may  be  seen  that  the  rate  coefficients  thus 
calculated  are  too  low  to  account  for  the  experimental  data  below 
about  1400  K.  If  a  higher  barrier  for  the  reverse  ruction  were 
to  be  included,  the  calculated  rate  coefficients  would  decrease 
further,  reinforcing  this  conclusion.  We  thus  speculate  that  at 
least  below  1400  K  the  main  products  of  reaction  2  arc  adducts 
correlating  with  the  singlet  and  triplet  potential  energy  surfaces 
of  AlCl*^  with  CO2'. 

The  AlO  +  SO2  Reaction.  The  abstraction  reaction 

AlCl  +  SO2  --  OAlCl  +  SO  (3a) 

has  an  endothermicity^^  of  119  kJ  moh^  and  is  spin-forbidden. 

For  temperatures  between  800  and  1100  K  the  rate  coefficient 
for  abstraction  is  much  lower  than  the  upper  limit  determination 
(cf.  eq  8)  because  of  the  large  endothermicity.  There  may  be 
channels  leading  to  the  formation  of  adducts.  However,  since  no 
reaction  could  be  detected,  the  adducts  must  have  low  binding 
energies  and  the  addition  channels  may  have  significant  activation 
barriers. 

The  AIG  +  O2  Reaction.  In  the  study ^  of  the  rate  coefficients 
of  the  AlCl  +  O2  reaction  from  490  to  1750  K,  it  was  speculated 
that  the  dominant  low  temperature  channel  could  be 

AlCl  +  O2  —  AlO.  +  Cl  (16a) 

while  at  higher  temperatures  the  channel 

AlCl  +  02^  OAlCl  +  0  (16b) 

was  considered  a  possibility.  Thermochemical  data  available  at 
that  time  indicated  that  the  thermochemistry  of  channel  1 6a  was 
uncertain,  while  channel  16b  was  thought  to  be  exothermic.^ 
Recent  high-level  ab  initio  computations'^  on  AIO2  now  indicate 
that  channel  16a  is  87  kJ  moh^  endothermic,  confirming  the 
XHfiAlOj)  assumed  in  the  JANAF  tables.^  That  theoretical 
work^^  also  suggests  that  channel  16b  is  66  kJ  mol"^  endothermic. 
These  data  and  the  use  of  MTST  allow  a  reevaluation  of  the 
measurements  of  the  AlCl  +  O2  reaction. 

The  parameters  for  the  molecular  model  as  well  as  the  rotational 
constants  and  vibrational  frequencies  of  the  transition  state  for 
channel  I6a  are  summarized  in  Table  V.  The  geometry  and 
force  constants  have  been  obtained  based  on  the  assumption  that 
the  transition  state  resembles  the  products  OAIO  +  Cl,  with  the 
halogen  atom  attacking  the  A1  atom  by  a  side-on  approach.  We 
assume  the  AlO  distance  and  force  constants  Fr(AlO)  and  Fr 
(OAIO)  for  the  transition  state  are  equal  to  those  of  AIO2.®  The 
remaining  parameters  have  been  obtained  similarly  to  those  for 
the  AlCl  +  N2O,  CO2  reactions.  If  the  barrier  height  for  the 
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TABLE  V:  Transition-Sute  Parameters  for  AlO^CP 
geometrical  parameter^  force  constant^ 


r(AlCl) 

2.34 

F^AICI) 

-2.0 

r(AlO) 

1.69 

F^AiO) 

4.2 

5(OA10) 

160 

F^(OAIO) 

0.24 

6(0A1C1) 

100 

F,(OAiCl) 

0.2 

r(A102Cl) 

0 

Ft(A102C1) 

0.02 

rotational  constants  {B1B2B3):  172  X  1()“5  cm"^ 
vibrational  frequencies  o/:  119,  176,  200, 678,  981  cm~‘ 

“  We  assume  a  C21,  arrangement  for  the  atoms.  The  parameters  of  the 
transition  state  are  based  on  the  products  of  the  endothermic  reaction 
AlCl  -I-  O2  OAIO  +  Cl.  *  Units  of  r  arc  10-^°  m;  units  of  d  and  t  arc 
deg.  ^  Units  of  Fr  are  10^  N  m“^;  units  of  F$  and  Fr  are  10"^*  N  m  rad*^. 


JTABLE  VI:  Transition-State  Parameters  for  CIAIQO^ 

geometrical  parameter"’*  force  constant^ 


KClAl) 

2.08 

F^ClAl) 

3.2 

r(AlO) 

1.63 

F.(A10) 

4.2 

7(00) 

1.72 

FAOO) 

-2.0 

6(C1A10) 

160 

F#(C1A10) 

0.3 

6(AlOO) 

90 

F,(A100) 

0.6 

t(C1A100) 

0 

F,(C1A100) 

0.02 

rotational  constants  (BiBzBi):  3.27  X  10"^  cm‘^ 
vibrational  frequencies  i?,:  125, 131, 245,  522,  907  cm"* 

®  We  assume  a  trans  arrangement  for  the  atoms.  The  parameters  of 
the  transition  state  arc  based  on  the  products  of  the  endothermic  reaction 
AlCl  -}-  O2  O  +  OAlCl.  *  Units  of  r  arc  10“*®  m;  units  of  8  and  r  arc 
deg.  ^  Units  of  F,  arc  10^  N  m"*;  units  of  Fj  and  Fr  arc  10"**  N  m  rad"^. 

T  (K) 

1800  1000  700  500 


Figure  3.  Arrhenius  plot  of  experimental  rate  coefficients  for  the 
previously*  studied  AlCl  +  O2  (open  circles)  reaction.  Rate  coefficient 
expressions  for  the  AlCl  +  02“^  OAlCl  +  O  (full  line)  and  AlCl  +  O2 
—  AIO2  +  Cl  (dashed  line)  channels  derived  using  the  MTST  method 
are  shown  for  comparison  with  experiments. 

reverse  reaction  is  neglected,  MTST  calculations  give 

*^MTST  =  2.75  X  10-'*  (r/K)'  *‘  X 

exp(-10460  K/7^  cm*  molecule"' s"'  (17) 

As  can  be  seen  in  Figure  3  the  MTST  calculations  given  by  the 
dashed  line  indicate  that  channel  1 6a  is  insignificant  at  all 
temperatures. 

The  MTST  parameters  for  channel  16b  are  summarized  in 
Table  VI  and  have  been  obtained  similarly  to  those  for  the  AlCl 
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+  CO2  reaction.  If  we  assume  a  negligible  barrier  height  for  the 
reverse  reaction,  the  MTST  expression  obtained  is 

exp(-7938  K/T)  cm^  molecule'^  s'*  (18) 

In  Figure  3  the  experimental  rate  coefficient  measurements  are 
shown  along  with  those  from  eq  18.  This  comparison  suggests 
that  channel  16b  can  be  of  significance  above  about  1500  K  in 
agreement  with  the  earlier  speculation.^ 

Recent  ab  initio  calculations^^  have  indicated  the  existence  of 
a  stable  adduct  AICIO2,  which  correlates  with  the  reactants  AlCl 
+  O2.  Hence,  the  dominant  reaction  pathway  at  temperatures 
below  1 7 50  K  is  likely  to  be  addition  with  the  abstraction  reaction 
1 6b  becoming  of  greater  importance  with  increasing  temperature. 
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Rate  coefficients  for  the  title  reactions  have  been  measured  in  a  high-temperature  fast-flow  reactor  (HTFFR). 
Weighted  fits  give  for  BCl  +  SO2  Jk(470-1690  K)  =  7.7  X  10-“(r/K)3O  exp(-1680  K/D  cm^  molecule"'  r' 
and  for  BCl  +  N2O  A:(670-1000  K)  =  5.1  X  lO-'^  exp(-6340  K/D  cm^  molecule"'  s"'.  2a  accuracy  limits 
are  about  ±25%.  Model-based  transition-state  theory  (MIST)  is  used  to  calculate  rate  coefficients  for  the 
0-atom  abstraction  channels  of  these  reactions,  as  well  as  for  the  previously  measured  BCl  reactioiis  vrith  O2 
and  CO2.  Approximate  agreement  between  theory  and  experiment  is  obtained  for  the  four  reactions,  indicating 
abstraction  to  be  dominant  for  each  reaction.  The  reactions  are  compared  to  the  corresponding  AlCl  reactions. 


Introdactioii 

In  previous  work,  the  high- temperature  fast-flow  reactor 
(HTFFR)  technique  has  been  used  to  make  rate  coefficient 
measurements  on  the  reactions  BCl  +  O2,  in  the  540- 1670  K 
range,*  and  BCl  +  CO2,  over  the  770-1830  K  range.^  Here  we 
extend  this  scries  by  reporting  measurements  on  the  reactions 

BCl  +  SO2  products  ( 1 ) 

BCl  +  NjO  —  products  (2) 

in  the  respective  470-1 690  K  and  670-1000  K  temperature  ranges. 

In  another  recent  paper,  we  reported  on  the  AlCl  reactions 
with  these  same  four  oxygen-containing  oxidants  and  developed 
a  variant  of  classical  transition-state  theory  (CTST).^  This  model- 
based  transition-state  theory  (MTST)  is  applicable  to  reactions 
with  early  activation  barriers  and  is  particularly  useful  when  the 
input  data  for  the  CTST  methods  originated  by  Benson^  are  not 
available.  MTST  was  used  in  that  work^  to  determine  the  AT^ 
part  of  the  rate  coefficient  k{T)  =  AT^  expf-E/  RT)  expressions 
for  abstraction  reactions.  Here  we  continue  the  comparison 
between  experiment  and  MTST  calculations  for  reactions  of  BCl, 
isoelectronic  with  AlCl. 

Experimental  Section 

The  HTFFR  technique,  as  used  for  BCl  studies,  has  been 
described  in  detail  elsewhere,^^  Briefly,  BCl  was  generated  by 
passing  20-40  ppm  BC^/Ar  through  a  microwave  discharge  in 
a  Pyrex  tube  (1.2-cm  i.d.)  and  fed  into  the  reaction  tube.  The 
reactor  consisted  of  a  vertically  positioned  mullite  tube  (2.2-cm 
i.d.)  sealed  inside  a  water-cooled,  steel  vacuum  housing.  Two 
parallel  series  of  SiC  heating  elements  were  used  to  maintain 
uniform  temperatures  in  the  reaction  tube.  Efficient  heating  of 
the  bath  gas  was  achieved  by  fixing  the  smaller  mullite  tube 
(1.8-cm  o.d.)  in  the  reaction  tube  upstream  from  the  reaction 
zone.  The  oxidant  gas,  SO2  or  N 2O  mixed  with  N 2,  was  introduced 
into  the  reaction  tube  through  a  movable  alumina  ring  inlet  with 
eight  circular  openings  (0.1 -cm  diameter)  through  which  the 
oxidant  gas  mixtures  were  injected  into  the  BCl/ Ar  stream.  Gas 
flow  through  the  inlet  was  maintained  at  a  constant  molar  rate 
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for  rate  coefficient  determinations,  each  of  which  consisted  of 
five  or  six  measurements  at  different  oxidant  concentrations.  At 
10  or  20  cm  downstream  from  the  ring,  relative  concentrations 
of  BCl  were  monitored  by  laser-induced  fluorescence  (LIF)  using 
four  different  transitions  of  the  (A*n-X*  2+)  band,  namely,  ( 1 ,0)Q 
at  266.0  nm,  (0,0)Q  at  272.0  nm,  (0,0)P  at  272. 1 6  nm,  and  ( 1 , 1  )P 
at  272.2 1  nm.^  Laser  pulses  at  these  wavelengths  were  generated 
byaLambdaPhysikEMG  101  excimer  laser/ FL  2002  dye  laser/ 
KDP  doubling  crystal  combination.  The  pulses  were  directed 
through  two  4-inm-diamcter  holes  in  the  mullite  reaction  tube. 
The  BCl  fluorescence,  which  passed  through  a  third  hole  in  the 
reaction  tube  positioned  normal  to  the  laser  pulse  path,  was 
recorded  with  an  EMI  9813  QA  photomultiplier  tube  (PMT)/ 
LeCroy  VVIOOB  wide-band  pulse  amplifier/Data  Precision 
Analogic  6000/620  100-MHz  transient  digitizer  combination. 
A  270-nm  interference  filter  (24-nm  fwhm)  was  positioned  in 
front  of  the  PMT  to  remove  background  radiation. 

Rate  coefficient  measurements  were  made  in  the  stationary 
inlet  mode"^  with  reaction  zone  lengths  of  10  or  20  cm,  under 
pseudo- first-order  conditions  with  BCl  as  the  minor  reactant. 
Plots  of  In  [BCllrelauve  VS  [SO2]  Or  [N2O],  in  which  [BCl]rd»tivt 
is  the  lime-integrated  LIF  signal  averaged  over  200  laser  pulses, 
were  linear  with  slopes  proportional  to  the  rate  coefficients.  For 
each  individual  measurement,  ki  and  <Tki  were  determined  by  a 
weighted  linear  least-squares  fit  of  the  data.”^’* 

The  gases  used  were  Ar  (99.998%,  Linde)  and  N2  (99.998%, 
Airco)  obtained  from  the  liquid,  BCI3  (99.9%,  Matheson),  N2O 
(99.99%,  Matheson),  N2O  (99.0%,  Linde),  SO2  (99.98%,  Mathe¬ 
son),  SO2  (99.98%,  Linde),  and  O2  (5.0%  in  Ar,  certified;  Scott). 

Results 

The  measured  rate  coefficients  fc,  of  reactions  1  and  2  at  the 
various  temperatures  used  arc  shown  in  Figure  1.  They  arc 
summarized  in  Tables  I  and  II  along  with  the  observed  reaction 
zone  length,  total  pressure  P,  total  concentration  [M],  maximum 
oxidant  concentration  [S02]™or  [N20]n«x,  fluorescence  intensity 
F  (a  measure  of  [BCl]  at  the  start  of  the  reaction  zone),  and 
average  velocity  v.  Variation  in  the  experimental  conditions  was 
used  to  check  the  accuracy  of  the  ki  measurements,  as  secondary 
reactions  of  potential  other  species  produced  in  the  Ar  +  BCl 
discharge,  i.e.,  B,  BCl,  BCI2,  Cl,  and  CI2,  could  have  produced 
BCl  in  the  reaction  zone.  If  this  had  occurred,  k,  measurements 
would  have  shown  dependences  on  reaction  zone  length,  D,  and 
P.  Our  reported  fc,  values  (Tables  I  and  II)  may  be  seen  to  be 
independent  of  these  quantities  by  examining  plots  of  [k{7)  - 
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Figure  1.  Arrhenius  plots  of  experimental  rate  coefficients  for  the  BCI 
+  SO2  (open  circles)  and  BCI  +  N2O  (open  squares)  reactions.  The  full 
lines  represent  the  experimental  rate  coefficient  expressions. 

ki]/k(T),  where  k{T)  are  the  fitted  rate  coefficient  expressions 
given  below,  vs  the  experimental  parameters;  for  each  case,  the 
it,  data  are  randomly  scattered  about  the  zero  line  of  the  ordinate. 

BO  +  SO2  Reaction.  The  85  rate  coefficient  measurements 
of  this  reaction  were  obtained  over  the  temperature  range  470-' 
1690  K  (Table  I).  In  this  study,  LIF  radiation  was  observed  in 
the  absence  of  BCI  at  temperatures  above  about  1 300  K.  For 
each  measurement  with  BCI,  the  intensity  of  this  radiation  was 
subtracted  from  the  signal.  The  interference  was  observed  to 
increase  as  the  reaction  tube  temperature  was  raised  and  as  the 
[SO2]  increased.  Since  SO2  as  well  as  SO,  generated  from  thermal 
dissociation  of  SO2,  absorb  in  the  260-“280-nm  region,^  both  could 
have  contributed  to  this  background  radiation.  There  was  also 
an  a  priori  possibility  that  SO,  produced  by  reaction  1,  could 
have  interfered  with  the  rate  coefficient  measurements.  This 
would  have  resulted  in  decreased  ki  values.  However,  since  the 
measured  rate  coefficients  are  independent  of  the  [SO2]  ranges 
and  the  initial  [BCI]  concentrations,  cf.  Table  I,  and  since  plots 
of  In  [BCljrriidveVS  [SO2]  are  linear,  wc  concludc  that  fluorescencc 
from  SO  did  not  significantly  affect  the  measurements.  Fur¬ 
thermore,  BCI  was  monitored  at  four  different  wavelengths,  cf. 
Table  I,  all  of  which  gave  similar  results. 

The  ki  measurements  were  fitted  to  the  expression  k(T)  - 
AT^  exp(-£//?T)  using  a  three-parameter  nonlinear  regression 
analysis.  Tlie  value  of  n,  determined  from  the  fit,  was  rounded 
to  one  decimal  place,  and  a  two-parameter  nonlinear  regression 
analysis  was  done.  The  resulting  expression  is 

A:, (470-1690  K)  =  7.7  X  lO-^^Cr/K)^®  X 

exp(-1680  K/71  cm^  molecule”'  s”'  (3) 
with  variance  and  covariance  matrix  elements^  ^ 

=  0.00142,  =  1164  K^  =  1.193  K 

(4) 

This  method  is  advantageous,  because  it  prevents  round-off  errors 
in  the  value  of  n  in  the  fmal  fit  and  simplifies  the  error  analysis. 
Combining  the  covariance  matrix  elements  by  the  propagation 
of  errors  technique  yields  ^2<Xk  precision  limits  of  8%  at  470  K, 


3%  at  1080  K,  and  4%  at  1690  K.  Allowing  10%  uncertainty 
associated  with  the  How  profile  factor^  and  a  conservative  20% 
for  further  systematic  errors  yields  ±2<Jk  confidence  limits  of 
24%  at  470  K  and  23%  at  1080  and  1690  K. 

BCI  +  N2O  Reaction.  Preliminary  measurements  showed  a 
dependence  of  kz  on  the  average  flow  velocity,  indicating  that 
continued  BCI  production  in  the  reaction  zone  was  occurring. 
Addition  of  50-500  ppm  of  O2  just  downstream  from  the 
microwave  discharge  removed  this  dependence,  and  this  procedure 
was  therefore  adopted.  In  Table  II,  the  rate  coefficient  mea¬ 
surements  on  reaction  2  over  the  temperature  range  670-1000 
K  arc  summarized.  Above  about  ICiOO  K  the  measured  rate 
coefficients  were  observed  to  be  dependent  on  pressure,  suggesting 
that  decomposition  of  N2O  is  interfering  with  the  measurements. 
Hence  these  are  not  included  in  the  table.  The  kz  measurements 
were  fitted  to  the  Arrhenius  expression  k{T)  =  A  cxp{-E/RT) 
to  give 

A:j(67(>-1000  K)  =  5.1  X  10”'^  X 

exp(-6340  K/T)  cm^  molecule”'  s”'  (5) 
with  matrix  elements 

=  0.0764,  (7£^  =  56400  K^  (r^/>l  =  65.1K  (6) 

Combining  these  matrix  elements  by  the  propagation  of  errors 
technique  yields  ±2(Tk  precision  limits  of  18%  at  670  K,  7%  at 
830  K,  and  10%  at  1000  K  and  ±2(Tk  confidence  limits  of  28% 
at  670  K,  24%  at  830  K,  and  25%  at  1000  K. 

Model-Based  Transition-State  Theory  (MTST) 

Recently, 3  we  developed  a  method  for  applying  transition- 
state  theory  to  reactions  of  metal-containing  species.  As  in  CTST, 
the  rate  coefficient  is  given  by^^ 

exp(-A£V^r»  (7) 

where  A£*  is  the  difference  in  zero-point  levels  of  the  transition 
state  and  the  reactants,  q*  is  the  partition  function  for  the  transition 
state,  and  and  are  partition  functions  for  the  reactants.  In 
the  MTST  method,  we  use  molecular  models  to  estimate  the 
vibrational  frequencies  and  rotational  constants  that  are  needed 
to  evaluate  q*.  For  the  molecular  model,  we  assume  that  harmonic 
restoring  forces  oppose  changes  in  bond  lengths,  angles  between 
valence  bonds,  and  dihedral  angles.  The  potential  energy  V  is 
thus  given  by^^ 

2K  =  («) 

1  J  * 

where  £r,  £«,  and  Fr  are  force  constants  for  changes  in  internuclear 
distances  r,  valence  angles  6,  and  dihedral  angles  r  from  their 
respective  equilibrium  values.  As  described  previously,^  the 
MTST  approach  is  suitable  for  reactions  with  transition  states 
which  arise  from  an  early  crossing  of  neutral  and  ionic  potential 
surfaces,  particularly  for  many  oxidation  reactions  of  metal- 
containing  species.  For  such  transition  states,  one  can  make 
reasonably  good  estimates  of  the  geometrical  parameters  and 
force  constants  based  on  Hammond’s  postulate.^^  For  example, 
for  an  exothermic  reaction  it  can  be  assumed  that  the  internuclear 
distances,  valence  angles,  and  force  constants  of  the  transition 
states  are  closely  related  to  those  properties  of  the  reactants.  The 
parameters  required  by  the  treatment  that  cannot  be  obtained 
from  properties  of  the  reactants  are  estimated,^  such  that  the 
calculated  vibrational  frequencies  of  transition  states  are  in  good 
agreement  with  ab  initio  computations* on  similar  reactions. 
We  have  shown**  that  for  transition  states  of  eight  reactions,  one 
of  the  reactants  being  a  B-  or  Al-containing  species  in  each  case, 
calculated  vibrational  frequencies  show  an  average  difTcrence 
with  ab  initio  results  of  less  than  about  100  cm**.  With  such 
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TABLE  I:  g-'nmary  of  Rate  Coefficient  Measurenwnts  for  BCl  +  SO^* 


rcation 

zone 

length 

(cm) 

p 

(mbar) 

cm"3) 

(10^5 

cm-3) 

c 

(ms-') 

T 

(K) 

k  ±<Tk 

(cm^  molecule-'  r') 

20 

40.5 

2.2 

0.6 

25" 

35 

1308 

4.35  ±  0.33  (-13)^ 

10 

40.6 

2.2 

1.1 

19" 

36 

1316 

4.80  ±0.35  (-13) 

10 

64.1 

3.6 

1.2 

55" 

34 

1295 

4.91  ±0.37  (-13) 

20 

64.0 

3.6 

0.6 

49" 

34 

1296 

3.93  ±0.32  (-13) 

20 

27.9 

1.2 

0.5 

27" 

91 

1622 

1.33  ±0.12  (-12) 

10 

27.9 

1.2 

1.0 

38" 

92 

1637 

1.65  ±0.13  (-12) 

20 

48.7 

2.2 

0.4 

21" 

53 

1626 

1.38  ±0.13  (-12) 

10 

48.7 

2.2 

0.8 

40" 

52 

1623 

1.40  ±0.11  (-12) 

10 

31.7 

3.2 

5.3 

33" 

18 

723 

3.62  ±0.36  (-14) 

20 

31.6 

3.3 

5-4 

15" 

18 

702 

2.43  ±0.23  (-14) 

20 

82.4 

9.6 

8.0 

27" 

12 

621 

1.13  ±0.08  (-14) 

20 

31.2 

1.4 

1.0 

37" 

101 

1628 

1.11  ±0.09  (-12) 

20 

31.6 

1.4 

1.0 

27" 

101 

1651 

1.15  ±0.09  (-12) 

10 

53.9 

2.4 

0.8 

24" 

66 

1650 

1.46  ±0.16  (-12) 

10 

62.0 

2.7 

0.7 

34 

1685 

1.24  ±0.11  (-12) 

10 

29.3 

1.3 

1.3 

17/ 

72 

1691 

1.45  ±0.13  (-12) 

20 

29.5 

1.3 

0.7 

ly 

71 

1691 

1.56  ±0.12  (-12) 

20 

33.5 

2.7 

3.4 

29 

28 

903 

8.02  ±  0.72  (-14) 

20 

67.6 

5.5 

1.1 

17/ 

18 

895 

7.69  ±0.58  (-14) 

20 

22,5 

1.8 

2.7 

AAf 

55 

904 

1.08  ±0.11  (-13) 

20 

23.6 

1.8 

2.9 

24 

33 

953 

'  03  ±  0.09  (-13) 

10 

46.1 

3.5 

2.9 

12 

17 

962 

0  68  ±  0.78  (-14) 

10 

28.6 

2.1 

2.6 

21 

37 

977 

1.35  ±0.14  (-13) 

10 

44.9 

3.3 

3.2 

20 

30 

975 

1.19  ±0.10  (-13) 

10 

78.9 

5.7 

2.8 

14 

17 

999 

1.19  ±0.10  (-13) 

10 

27.6 

1.5 

1.5 

14 

66 

1377 

5.99  ±0.91  (-13) 

10 

85.2 

4.5 

1.0 

16 

29 

1368 

4.62  ±0.48  (-13) 

10 

54.2 

2.9 

1.2 

24 

40 

1375 

4.91  ±0.49  (-13) 

10 

54.2 

2-9 

1.2 

4" 

40 

1375 

5.26  ±0.40  (-13) 

20 

25.7 

1.4 

0.6 

9" 

82 

1329 

5.23  ±0.69  (-13) 

20 

25,7 

1.4 

1.2 

39 

82 

1329 

5.53  ±0.54  (-13) 

20 

59.1 

3.3 

0.5 

76 

42 

1311 

4.89  ±0.39  (-13) 

20 

33.5 

1.8 

0-6 

88 

50 

1346 

5.37  ±0.45  (-13) 

20 

20.1 

1.0 

1.0 

50 

98 

1472 

9.18  ±0.82  (-13) 

20 

62.2 

3.2 

0.6 

68 

49 

1430 

7.36  ±0.59  (-13) 

20 

39.5 

1.9 

0.8 

44 

60 

1507 

9.35  ±0.81  (-13) 

20 

34.5 

1.6 

0.5 

31 

59 

1554 

1.06  ±  0.07  (-12) 

10 

34.6 

1.6 

1.6 

34 

60 

1570 

8.37  ±0.60  (-13) 

10 

94.6 

4.5 

1.2 

63 

35 

1535 

7.33  ±0.53  (-13) 

10 

58.2 

2.6 

0.9 

20 

44 

1602 

1.38  ±0.16  (-12) 

10 

79.9 

3.8 

0,7 

15 

31 

1536 

1.09  ±0.11  (-12) 

10 

19-7 

1.4 

2.5 

31 

38 

1003 

1.69  ±0.19  (-13) 

10 

41.1 

2.9 

1.8 

38 

27 

1041 

2.10  ±0.18  (-13) 

rcation 

zone 

length 

(cm) 

P 

(mbar) 

[Ml 

(10” 

cm"S) 

[S02]tm. 

(10'5 

cm“^) 

(m  r') 

r 

(K) 

k±9t 

(cm^  molecule-'  s“') 

10 

55.3 

3.9 

3.3 

48 

29 

1027 

1.89  ±0.12  (-13) 

10 

66.3 

4.5 

1.6 

35 

19 

1073 

2.20  ±0.22  (-13) 

20 

46.7 

3.1 

l.l 

34 

27 

1080 

2.00  ±0.18  (-13) 

20 

25.9 

1,6 

0.8 

32 

35 

1138 

2.65  ±0.25  (-13) 

20 

30.5 

2.0 

1.7 

53 

57 

1102 

2.39  ±0.17  (-13) 

20 

39.6 

2.5 

I-i 

32 

35 

1126 

2.46  ±0.26  (-13) 

10 

39.7 

2.4 

2.6 

56 

37 

1187 

3.61  ±0.27  (-13) 

10 

89.3 

5.5 

1.4 

45 

21 

1171 

3.02  ±0.18  (-13) 

10 

50.3 

3.0 

1.5 

43 

21 

1212 

3.01  ±0.22  (-13) 

10 

32.1 

1.8 

2.3 

66 

42 

1258 

4.49  ±  0.34  (-13) 

20 

32.0 

1.8 

1.2 

33 

43 

1277 

3.77  ±  0.30  (-13) 

20 

44.1 

2.6 

1,7 

52 

46 

1226 

3.45  ±  0.22  (-13) 

20 

19.9 

1.8 

3,2 

37 

30 

802 

4.99  ±0.51  (-14) 

20 

60.7 

5,8 

2.6 

17 

16 

762 

4.20  ±0.30  (-14) 

20 

39.7 

3.9 

5.0 

29 

19 

746 

3.33  ±0.25  (-14) 

20 

49.0 

5.1 

4.4 

13 

11 

696 

2.32  ±0.27  (-14) 

10 

49.1 

5.1 

8.7 

15 

11 

692 

2.48  ±0.16  (-14) 

10 

83.2 

9.6 

6.2 

13 

8 

627 

1.41  ±0.11  (-14) 

20 

28.5 

4.0 

9.6 

11 

10 

523 

5.40  ±0.50  (-15) 

20 

79.8 

12.3 

15.5 

12 

6 

471 

2.12  ±0.19  (-15) 

20 

63.3 

9.6 

15,9 

19 

6 

476 

2.77  ±0.17  (-15) 

20 

63.3 

9.1 

15.2 

18 

6 

503 

2.87  ±0.26  (-15) 

20 

57.7 

8.5 

10.5 

53 

9 

490 

2.31  ±0J2(-1S) 

20 

85.4 

12.2 

15.0 

7 

6 

508 

4.05  ±  0.53  (-15) 

20 

63,9 

8.8 

14.0 

24 

7 

527 

4.65  ±0.46  (-15) 

20 

35.4 

4.7 

11.1 

12 

9 

548 

7.37  ±0.62  (-15) 

20 

56.1 

7.2 

12.7 

17 

8 

563 

6.48  ±0.39  (-15) 

20 

48.4 

6.1 

10.7 

25 

9 

575 

8.00  ±  0.62  (-15) 

20 

54.7 

7.2 

8.4 

38 

11 

552 

4.94  ±0.60  (-15) 

20 

47.1 

5.7 

8.2 

36 

12 

600 

1.16  ±0.09  (-14) 

20 

55.6 

6.9 

7.9 

31 

12 

583 

7.22  ±0.50  (-15) 

10 

85.3 

10.4 

11.9 

21 

8 

592 

1.05  ±0.13  (-14) 

10 

47.6 

5.2 

7.3 

55 

13 

668 

2.45  ±0.20  (-14) 

10 

77.1 

8.8 

8.7 

34 

11 

635 

1.79  ±0.23  (-14) 

10 

40,9 

4.5 

5.6 

93 

17 

657 

2.34  ±0.26  (-14) 

20 

20.0 

2.1 

5.3 

84 

18 

674 

1.73  ±0.19  (-14) 

20 

24.8 

2.6 

6.1 

21 

16 

678 

1.98  ±0.16  (-14) 

20 

20.6 

1.7 

3.4 

37 

28 

897 

8.38  ±0.72  (-14) 

20 

51.1 

4.3 

4.3 

51 

22 

854 

5.48  ±0.33  '  4) 

20 

39.5 

3.3 

4.7 

32 

20 

872 

6.12  ±0.48  (-14) 

10 

39.5 

3.3 

4.8 

46 

20 

859 

7.04  ±0.68  (-14) 

10 

76.2 

6.7 

3.3 

32 

15 

826 

6.93  ±0.55  (-14) 

•  The  measurements  are  reported  in  the  sequence  in  which  they  were  obtained.  ‘  ^  arbitra^  um«.  Unl^  othe,^  “Si  wS 

at  the  272.0-nm  transition.  BCl  was  monitored  at  the  272.16-nm  transition.  ^  Should  be  read  as  (4.35  ±  0.33)  X  10-  cm  molecule-  s-  .  BCl  was 
monitored  at  the  266.0-nm  transition.  ^  BCl  was  monitored  at  the  272.21-nm  transition. 


estimated  geometrical  parameters  and  force  constants,  we  obtain 
the  vibrational  frequencies  by  a  normal  coordinate  analysis  of 
the  molecular  model,  which  procedure  has  been  described  in  detail 
by  Pauling  and  Wilson.*^  Finally,  rotational  constants  are 
calculated  from  the  eigenvalues  of  the  inertia  tensor.^* 

Discussion 

BQ  +  SO2  Reacdofl.  The  thermochemically  accessible 
channels  of  reaction  1 

BCl  +  SO2  --  OBCl  +  SO  (la) 

--BO2  +  SCI  (lb) 

are  245  and  26  kJ  moh*  exothermic,  respectively,  at  0  K.^^22 
Channel  lb  is  probably  inhibited  by  a  substantial  energy  barrier, 
as  all  of  the  reactant  bonds  must  be  broken  and  reformed.  In 
contrast,  we  expect  channel  la  to  have  a  slight  barrier  as  only 
the  O-SO  bond  is  broken  in  the  reaction.  In  channel  1  a,  formation 
of  SO(^2'),  the  ground  state,  is  forbidden  by  the  spin  rule;  however, 
formation  of  SO(^  is  both  exothermic  (by  169  kJ  mohO  and 
spin  allowed.  The  MTST  calculations  which  follow  are  based  on 
the  formation  of  this  state.  Table  III  summarizes  the  rotational 
constants  and  vibrational  frequencies  for  the  transition  state.  The 


first  two  terms  of  cq  7  are  computed  at  several  temperatures  in 
the  range  500-2000  K  and  fitted  to  the  expression  ^7^.  An 
estimate  of  the  barrier  height  is  made  by  equating  eqs  3  and  7 
at  735  K,  the  midpoint  of  the  range  investigated,  which  yields 
24.1  kJ  moh^  The  present  calculations  lead  to  the  expression 

*MTsr  =  519  X  10-'’(r/K)'  “  X 

exp(-2902  K/7)  cm^  molecule"*  s"‘  (9) 

Figure  2  compares  the  MTST  calculations  with  the  recommended 
/t,(D  expression,  and  they  are  in  good  agreement  over  the  470- 
1690  K  range.  This  agreement  suggests  that  channel  la,  for 
which  the  MTST  calculations  are  based,  is  the  dominant  pathway. 
In  a  related  paper,^  we  attempted  to  make  rate  coefficient 
measurements  on  the  homologous  AICl  +  SO2  reaction.  We 
failed  to  observe  any  significant  disappearance  of  AlCl  by  reaction 
with  SO2,  indicating  that  the  AlCl  reaction  has  a  much  lower 
rate  coefficient  than  the  BCl  SO2  reaction.  The  lower  rate 
coefficient  is  consistent  with  the  reaction  being  at  least  1 10  kJ 
mol”'  endothermic.*-^ 

BQ  +  N2O  Reaction.  The  abstraction  pathway  of  reaction  2 

BCl  +  NjO  —  OBCl  +  Nj  (2a) 

is  631  kJ  mol-'  exothermic  at  0  K.*-*  Similarly  to  our 
investigation  of  the  above  treatment  of  BCl  +  SO2,  we  investigate 


Kinetics  of  BCl  +  Oxidant  —  Product 
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TABLE  m  Summary  of  Rate  Coefficient  Measurements  for  BO  -t- 


reaction  zone 
length  (cm) 


20 

20 

20 

20 

20 

20 

20 

20 

20 

20 

20 

20 

20 

20 

20 

20 

20 

20 

20 

20 

20 

20 

20 

20 

10 

10 

10 


P  (mbar) 


42.3 

36.3 

36.3 

84.6 

84.8 

40.7 

78.8 

59.1 

99.2 

99.4 

40.5 
40.5 

77.2 

55.1 

56.8 

56.8 
51.7 

51.7 
42.0 
42.0 

78.1 

78.1 

55.7 

55.7 

53.7 
53.7 

47.1 


[M] 

(10^'^  cm"^) 


4.5 

3.4 

3.5 

8.7 
8.9 
3.0 
6.1 

4.3 

7.3 
7.3 

3.8 

3.8 

7.3 
5.0 
5.0 
5.0 

4.3 

4.3 
3.2 
3.2 

5.9 
5.9 
5.9 
5.8 

4.5 

4.4 

3.6 


[NzOln^. 
(10^^  cm”^) 


V  (m  s"') 


23.2 
16.1 

16.4 

50.2 

51.4 

12.2 

27.8 

23.2 
16.1 
16.1 
20.7 
20.6 

39.3 

24.6 

28.7 

28.5 

19.8 
19.7 

17.6 

17.4 
22.2 
22.2 

40.2 

39.7 

27.2 
27.1 
23-0 


15 

58 

59 
27 
24 
35 
23 
42 

27 
26 
45 
35 

37 

42 
76 
79 

43 

38 
84 
72 
29 
33 

39 

28 
24 
20 
21 


13 

18 

18 

11 

11 

24 

20 

24 
18 
18 
21 
21 
11 
12 
20 
20 
15 

15 

25 
25 
13 
13 
11 
11 

16 
16 
25 


«  The  measurements  are  reported  in  the  sequence  in  wUch  were  obu^.  *  In  arbitrary  units, 
the  272.0-nm  transition.  ‘Should  be  read  as  (4.28  ±  0.40)  x  10-  cm  molecule 


r(K) 

k±  at 

(cm^  molecule”*  r*) 

674 

771 

4.28  ±  0.40  {-16Y 
1.94  ±0.34  (-15) 

756 

1.78  ±0.14  (-15) 

709 

5.50  ±  0.73  (-16) 

691 

4.64  ±0.31  (-16) 

982 

6.07  ±  0.59  (-15) 

939 

3.88  ±0.69  (-15) 

997 

7.14  ±0.64  (-15) 

989 

1.08  ±  0.07  (-14) 

989 

1.06  ±  0.09  (-14) 

768 

9.14  ±0.97  (-16) 

772 

1.80  ±0.42  (-15) 

770 

1.07  ±0.10  (-15) 

792 

1.48  ±0.12  (-15) 

822 

2.59  ±0.39  (-15) 

829 

2.89  ±0.28  (-15) 

872 

4.01  ±0.43  (-15) 

879 

4.42  ±0.40  (-15) 

945 

7.09  ±  0.60  (-15) 

956 

8.23  ±0.57  (-15) 

960 

5.56  ±  0.43  (-15) 

961 

5.74  ±0.39  (-15) 

687 

4.29  ±  1.04  (-16) 

696 

5.95  ±  1.06  (-16) 

873 

5.22  ±0.45  (-15) 

877 

4.27  ±  0.52  (-15) 

950 

5.45  ±  0.74  (-15) 

For  each  measurement,  BCl  was  monitored  at 


TABLE  ni:  Transition-State  Parameters  for  GBOSO* 


geometrical  parameter^ 


r(ClB) 

KBO)  1-94 

r(OS)  1-43 

KSO)  1-43 

fl(ClBO)  120‘ 

9(BOS)  120" 

9(OSO)  120* 

t(CIBOS)  0° 

t(BOSO)  90“ 


rotational  constants  (BiBjBj) 
vibrational  frequencies  Vi 


force  constant^  _ _ 

FACIB)  3.5 

fABO)  -2.0 

FAOS)  5.15 

FASO)  10.3 

fKClBO)  0.2 

F,(BOS)  0.6 

FAOSO)  1.68 

f,(ClBOS)  0.02 

f.(BOSO)  0.02 

7.95  X  10-^  cm"’ 

50, 65,  122,237,  520,839, 
896,  1300  cm-* 


"  We  assume  trans  arrangements  for  the  atoms  in  the  CIBOS  an 
BOSO  groups.  ‘  Units  of  r  are  lO"*®  m.  ‘  Uniu  of  F,  are  10^  N  m*  ;  units 
of  F,  and  F,  are  10-*«  N  m  rad'^. 


the  kinetics  of  this  reaction  path  using  MTST  with  the  parametere 

summarized  in  Table  IV.  The  first  two  terms  of  eq  7  arecomputed 
in  the  500-2000  K  range  and  fitted  to  the  expression  AT'.  Usmg 
eqs  5  and  7  at  the  mid  r*  point  802  K,  we  back  out  a  burner 
height  of  53.2  kJ  mol'*,  which  combines  with  the  calculated  A 
and  n  values  to  give 

=  7.98  X  lO-'^r/K)'®'  X 

exp(-6399  K/D  cm^  molecule"’  s’  (10) 

As  shown  in  Figure  2,  eq  10  is  roughly  in  agreement  with  the 
recommended  ki(,T)  expression  in  the  670-1000  K  temperature 
range.  We  obuined  similar  results  in  a  previous  study^  of  the 
homologous  AlCl  -l-  NjO  reaction.  The  two  reactions  are  cicely 
related  in  that  both  have  fairly  large  activation  barriers  and  both 
reactions  have  exothermic  channels  leading  to  O-atom  abstraction. 
BO  +  Oi  Reaction.  The  abstraction  reaction 

BCH-O.-^OBCH-O  (11) 

is  298  kJ  mol-*  exothermic  at  0  We  again 

transition-state  parameters,  which  are  summarized  in  Table  V, 
for  the  abstraction  channel.  By  equating  the  MTST  expression 


T  (K) 


re  2.  Arrhenius  plots  of  experimenul  rate  coemdents  for  ;^e  BCl 
.  \ _ 1  cnitarps^  rcactlOOS.  TuC  lull 


to  the  experimental  rate  coefficient  expression*  at  820  K,  we 
calculate  a  barrier  height  of  31.0  Id  mol-*,  which  yields 

W  =  5.i7xio-”(r/K)‘’’x 

exp(-3730  K/T)  cm’  molecule"’  s"’  (12) 

As  shown  in  Figure  3,  eq  12  is  in  good  agreement  with  the 
experimental  measurements.  In  contrast  to  the  present  study. 


Futcrko  et  al. 
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TABLE  IV:  Transition-SUte  Parameters  for  CIBONN*  TABLE  VI;  Transition-State  Parameters  for  GBOCO* 

geometrical  parameter^  force  constant*^  geometrical  parameter®**  force  constant^ 


/■(CIB) 

1.72 

FXCIB) 

3.5 

/IBO) 

1.94 

F^BO) 

-2.0 

KON) 

1.18 

F^ON) 

5.85 

r(NN) 

1.13 

F^NN) 

18.7 

«(CIBO) 

120® 

F^CIBO) 

0.2 

tf(BON) 

90® 

F^BON) 

0.6 

9(ONN) 

160® 

F^ONN) 

0.67 

t(CIBON) 

0® 

Fr(ClBON) 

0.02 

t(BONN) 

0® 

F.(BONN) 

0.02 

routional  constants  2.46  X  1 0-3  <^-3 

vibrational  frequencies  vi  58,  123,  170,  292,  596,  839, 

996, 2221  cm-^ 


®  We  assume  irans  arrangements  for  the  atoms  in  the  CIBON  and 
BONN  groups.  *  Units  of  r  arc  1 0"*°  m.  ^  Units  of  /v  are  1 CF  N  m“‘ ;  units 
of  Fs  and  Fr  are  10"^*  N  m  rad“^. 

TABLE  V:  Transition-State  Parameters  for  CIBOO^ 


geometrical  parameter^  force  constant^ 


r(ClB) 

1.72 

FXCIB) 

3.5 

riBO) 

1.94 

FXBO) 

-2.0 

riOO) 

1.21 

F.(00) 

5.9 

^(CIBO) 

120® 

FKCIBO) 

0.2 

^(BOO) 

120® 

F#(BOO) 

0.6 

t(CIBOO) 

0® 

Fr(ClBOO) 

0.02 

rotational  constants  {818182) 

1.00  X  10“2cm-^ 

vibrational  frequencies  V{  75, 131,  317,  841,  1116cm“^ 

®  We  assume  a  trans  arrangements  for  the  atoms  in  the  CIBOO.  *  Units 
of  r  arc  10”*°  m.  ^  Units  of  Fr  arc  10^  N  m-‘;  units  of  Fe  and  F^  arc  lO’^* 
N  m  rad‘2. 

r  {K) 


1700  1000  700  500 


0.6  1.1  1.6  2.1 


1000/  r  (K-^) 

Figure  3.  Arrhenius  plots  of  measurements  of  the  rate  coefficients  for 
the  BCl  +  0:  (open  squares)  and  BCl  +  CO2  (open  circles)  reactions 
from  previous  studies.  The  full  lines  represent  rate  coefficient 
expressions  for  the  0-atom  abstraction  channels,  as  calculated  by  the 
MTST  method. 

we  previously  found^  that  the  AlCl  +  O2  reaction  MTST 
calculations,  which  were  based  on  an  0-atom  abstraction 
mechanism,  were  inconsistent  with  experiments.  The  AlCI  +  O2 
reaction  differs  from  the  analogous  BCl  reaction  in  that  formation 
of  the  abstraction  products  OAICl  +  O  is  an  endothermic  process 
with  a  barrier^  of  at  least  66  kJ  mo\-K  On  the  basis  of  the  ab 
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1.72 
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3.5 

KBO) 

1.94 
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-2.0 
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1.16 

F^OC) 

8.0 

KCO) 

1.16 

PACO) 

16.0 

d(CIBO) 

120® 

FKCIBO) 

0.2 

i9(BOC) 

90® 

FAfiOC) 

0.6 

^(OCO) 

160® 

FAQCO) 

0.79 

r(ClBOC) 

0® 

F,(ClBOC) 

0.02 

r(BOCO) 

0® 

F,(BOCO) 

0.02 

rotational  constants  {81828^)  2.26  X  10-3  cin“3 

vibrational  frequencies  Wi  57, 122,  169,  303,  662,  842, 

1117, 2185  cm-* 

We  assume  trans  arrangements  for  the  atoms  in  the  CIBOC  and 
BOCOgroups.  *  Units  of  rare  10“*®  m.  ^  Units  of  F,  arc  lO^Nm”* ;  units 
of  F$  and  F,  are  10“**  N  m  rad-^. 

initio  study  by  Chen,  Hasc,  and  Schlegcl,*’  we  suggested  that  the 
main  channel  of  AlCl  +  O2  reaction  involves  the  formation  of  an 
adduct. 

BQ  +  CO2  Reaction.  The  O-atom  abstraction  reaction 

BCl  +  CO2  —  OBCl  +  CO  (13) 

is  266  kJ  mol”*  exothermic  at  0  The  calculated  transition- 

state  parameters  for  this  path  are  summarized  in  Table  VL  By 
equating  the  MTST  expression  to  the  experimental  rate  coefUdent 
expression^  at  1080  K,  we  calculate  a  barrier  height  of  73.2  kJ 
mol”*,  which  gives 

^MiST  =  1-23  X  10-^^(7yK)^®^  X 

exp(-^804  K/ 7)  cm^  molecule"^  s“^  (14) 

Figure  3  shows  a  comparison  of  eq  14  with  experimental  data 
from  the  previous  study  The  calculations  are  in  approximate 
agreement  with  experiment  over  the  770-1830  K  range.  A 
comparison  can  be  made  between  the  BCl  and  AlCl  reactions 
with  CO2.  As  is  the  case  for  the  BCl  and  AlCl  reactions  with 
O2,  the  BCl  +  CO2  reaction  has  an  exothermic  abstraction  channel 
while  the  channel  leading  to  formation  of  OAICl  +  CO  has  an 
enthalpy  barrier  of  at  least  98  kJ  mol“*.  Because  the  abstraction 
channel  for  the  AlCl  +  CO2  reaction  is  endothermic,  other 
channels  resulting  in  adduct  formation  can  dominate,  especially 
at  low  temperatures.^ 

Conclusions 

Kinetic  measurements  on  BCl  reactions  with  oxygen-containing 
oxidants  have  been  made  over  wide  temperature  ranges.  We 
used  a  TST- based  method  to  calculate  pre-exponentials  for  four 
BCl  reactions  and  obtained  good  agreement  with  experiments. 
We  showed  that  for  these  reactions  rate  coefficients  can  be 
reasonably  well  predicted  over  a  range  of  temperatures,  given  a 
rate  coefficient  measurement  at  a  single  temperature  to  obtain 
the  E  factor.  For  exothermic  processes,  i.e.,  ^e  N2O  reactions, 
AlCl  and  BCl  behave  similarly,  as  may  be  expected  for  cogeners. 
The  observed  differences  for  the  other  reactions  are  in  accord 
with  the  thermochemistry. 
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APPENDIX  E 


Activation  barriers  for  series  of  exothermic  homologous  reactions. 

IV.  Comparison  of  measurements  to  theory  for  reactions  of  atoms 
with  N2O 
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The  kinetics  of  the  AH-N20-*A10+N2  reaction  has  been  studied  in  a  high-temperature  fast-flow 
reactor  (HTpFR).  The  expression  k{520-\030  K)=4.6X10"’^  (T/K)^^  exp(-778  KJT) 
cm^  molecule"^  s"‘  has  been  derived  from  a  nonlinear  regression  analysis  of  the  measured  rate 
coefficients.  2a  precision  limits  are  about  ±5%  and  accuracy  limits  are  estimated  to  be  about 
±23%.  Combining  the  data  with  a  296  K  literature  measurement  yields  /:(296-l030 
K)=5.8X  10’'^  (T/K)‘  exp(-90.4  K/T)  cm^  molecule"'  s"',  with  2cr estimated  accuracy  limits  of 

about  ±27%.  The  semiempirical  approach  used  previously  to  predict  activation  barriers  for  5'  and 
5^  metal  atom  reactions  with  N2O  has  been  modified  to  allow  predictions  of  group  13  atom  reactions 
with  NoO.  The  activation  barriers  for  the  B,  Ga,  In,  and  Tl  atom  reactions  are  calculated  from  the 
experimental  activation  barrier  of  the  A1  +  N20  reaction.  ©  1995  American  Institute  of  Physics. 


I.  INTRODUCTION 

In  the  first  paper'  of  this  series  it  was  shown  that  the 
experimental  activation  barriers  E,  and  as  a  result  the  rate 
coefficients  of  reactions  of  ground  state  5'  and  s^  metal  at¬ 
oms  with  N2O,  are  closely  related  to  the  sums  of  the  s-p 
promotion  energies  and  ionization  potentials  of  the  atoms.  E 
was  defined  by 

k(T)  =  Ar  expi~E/RT)  (1) 

which  expression  gave  the  best  correlations  for  n=0.5.  A 
semiempirical  treatment  was  developed  to  calculate  E  for  all 
the  5'  and  5^  metal  atoms  reactions  with  N2O,  using  the 
experimentally  determ;  1  activation  barrier  of  the  CU  +  N2O 
reaction.  In  the  second  paper,'  the  treatment  was  refined  for 
5 '  metals  for  which  the  reaction  coordinate  had  to  be  taken 
as  a  variable  instead  of  a  constant.  In  the  third  paper,^  the 
treatment  for  5“  metals  was  further  modified  and  reactions 
with  O2  were  considered.  Here,  the  calculations  are  extended 
to  reactions  of  Z(5'p')  atoms  with  N2O,  where  the  bamer 
heights  for  the  group  are  determined  from  measurements  on 
the  AI  +  N2O  reaction,  reported  here. 

II.  EXPERIMENT 
A.  Technique 

The  measurements  on 

AH-N2O— AlO-t-N.  (2) 

for  which^  ^7^298  **  “345  kJ  mor',  were  performed  in  an 
HTFFR  (high-temperature  fast-flow  reactor).  This  facility^'^ 
and  the  operational  procedures^  have  been  described  in  de¬ 
tail.  Briefly,  a  vertical  mullite  (McDanel  MV- 30)  reaction 
tube  (60  cm  long,  2.2  cm  i.d.)  was  radiatively  heated  to  the 


‘^’Present  address:  Department  of  Chemical  Engineering.  Massachusetts  In¬ 
stitute  of  Technology,  Cambndge.  Mas.sachusetts  02139. 
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desired  temperature  by  columns  of  resistively  heated  SiC 
rods  inside  an  insulated,  water-cooled  vacuum  housing.  A1 
vapor  was  produced  by  resistively  heating  an  Al-wetted  W 
coil  in  a  flow  of  Ar  carrier  gas.  Mixtures  of  N2O  in  Ar,  at 
volume  flow  rates  of  between  2%  and  7%  of  the  main  Ar 
flow,  were  injected  into  the  Al/Ar  gas  stream  through  a  mov¬ 
able  inlet.  The  total  molar  flow  rate  of  gas  through  the  inlet 
was  kept  constant  during  a  given  rate  coefficient  measure¬ 
ment.  Further  downstream,  relative  Al  concentrations  were 
determined  by  laser-induced  fluorescence  (LIF)  or  atomic 
resonance  absorption  spectrometry  (ARAS).  ARAS  was  used 
as  a  check  on  the  LIF  results.  These  two  optical  measure¬ 
ment  techniques  have  an  important  difference.  LIF  measures 
the  concentration  of  Ai  at  the  axis  of  the  reactor,  whereas 
ARAS  measures  the  average  concentration  over  the  full 
width  of  the  reactor.  Comparison  thus  also  gives  an  indica¬ 
tion  of  the  effect,  if  any,  of  the  flow  profile.  For  the  LIF 
measurements  a  pulsed  Lambda  Physik  EMG  101  excimer/ 
FL2002  dye  laser  in  combination  with  a  KDP  doubling  crys¬ 
tal  was  used.  The  Al  {~S-^P°)  transition  at  265.2  nm  was 
pumped  and  the  fluorescence  intensity  was  measured  through 
a  262  nm  (26  nm  full  width  at  half-maximum)  filter  by  an 
EMI  9813  QA  photomultiplier  tube  (PMT),  connected  to  a 
Data  Precision  Analogic  6000/620  100  MHz  transient  digi¬ 
tizer.  For  the  ARAS  experiments  a  hollow  cathode  lamp  was 
used  to  generate  the  Al  line  radiation  at  309.3  nm  and  the 
transmitted  radiation  intensity  was  monitored  with  a  mono¬ 
chromator  in  combination  with  the  PMT  connected  to  a  Kei- 
thley  614  electrometer. 

Rate  coefficient  measurements  were  made  in  the  station¬ 
ary  inlet  mode^  at  observed  reaction  zone  lengths  of  10  or  20 
cm,  under  pseudo-first-order  conditions  with  Al  as  the  minor 
reactant.  [Al^iative  is  defined  as  [Al]/[All ;  where  [Al],  is  the 
concentration  of  Al  at  the  lowest  [N2O]  used.  The  [N2O]  was 
typically  varied  by  a  factor  of  5.  In  LIF,  [Al],.gimjyg^E/E, , 
where  F,  is  the  initial  fluorescence  and  F  the  fluorescence  in 
the  presence  of  N2O.  In  ARAS,  the  Lambert-Beer  law  gives 
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[AlLianv/'^in(/o//).  vvhere  /„  is  the  light  intensitv  passing 
through  the  reactor  without  Al  present  and  /  is  the  intensity 
with  A!  atoms.  Rate  coethcients,  k,,  were  obtained  trom  the 
slope  of  linear  plots  ot  ln[AllrL.};,i,vc  ^'s  [N2O].  A  weighted 
linear  regression  was  used  to  determine  k^  and  tor  each 
measurement/ 

The  gases  used  were  99.998%  Ar  from  the  liquid  Ar 
(Linde),  and  0.53%  N.O  (99.99%)  in  Ar  (99.999%)  from 
Matheson,  the  calibration  ot  which  was  checked  against  pure 

B.  Results  and  discussion 

The  measured  rate  coethcients  k,  and  the  experimental 
conditions  under  which  they  were  obtained  are  summarized 
in  Table  I.  The  parameters  varied  include:  total  pressure  P, 
corresponding  to  total  concentration  [M ],  initial  fluorescence 
intensity  F  or  initial  percentage  absorption  (both  a  measure 
of  initial  [Al]),  average  gas  velocity  J,  and  reaction  zone 
length.  Our  reported  k,  values  were  determined  to  be  inde¬ 
pendent  of  these  parameters  by  examining  plots  ot 
[k{T)~k^\fk(T).  where  A;(T)  is  obtained  from  Eq.  (3)  below 
vs  the  particular  parameter.  For  each  parameter  the  k^  data 
were  found  to  be  randomly  scattered  about  the  zero  line  ot 
the  ordinate.  Figure  I  shows  the  /c,  data  taken.  The  range  of 
temperatures  covered  is  520-1030  K.  The  lower  limit  was 
determined  by  the  heating  of  the  Ar  bath  gas  by  the  Al 
source.  The  upper  limit  was  set  by  the  onset  of  potential 
dissociation  of 

A  weighted  nonlinear  ht”  ot  the  data  to  Eq.  (1)  with 
n=0.5  yields 

A:(520-I030  K) 

=  4.6X  10-'"(7/K)'"-' 

Xexp( -6.469//?7')cm^  molecule”*  s  (3) 

where  E  is  expressed  in  kJ  mol  *.  The  covariance  matrix 
elements’^  are  calculated  to  be  =5.02x10  A  , 
cx^^=4.84XlO”'*  AE,  and  fTi  =  4.83x  lO”'  E\  The.se  yield 
±2(Tj^  precision  limits  varying  from  ±2%  to  ±7%.  Fitting 
only  the  LIF  results  w'ould  have  led  to  a  k{  7)  expres.sion 
2%-ll%  higher  than  Eq.  (3),  while  considering  only  the 
ARAS  data  yields  an  expression  10%-13%  low'er.  The  later 
set  consists  of  less  data  and  is  therefore  less  precise.  None¬ 
theless,  this  suggests  a  small  systematic  difference,  well 
within  the  ±20%  systematic  errors  and  ±10%  flow  profile 
factor  uncertainty  always  allowed  for  in  this  type  ot 
work.^’^’*^*  '^  Using  these  factt)rs  leads  to  accuracy  intervals 
varying  from  ±22%  to  in 24%.  The  ht  is  shown  in  Fig.  1. 

The  only  previously  reported  kinetic  measurements  on 
the  Al±NiO  reaction  resulted  in  a  k{296  K)  of 

(ll±l)xl0~'*-  cm' molecule  s  '.  It  was  obtained  by 
Mitchell  and  co-workers  using  their  static  cell  laser- 
photolysis,  laser-fluorescence  technique.'^  Their  uncertainty 
value  pertains  to  precision  only.  Allowing  tor  possible  sys¬ 
tematic  errors  of  similar  magnitude  as  in  our  work,  leads  to 
(1I±2  to  3)X10"'“  cm' molecule' ' s”'.  Extrapoladon  of 
Eq.  (3)  to  296  K  yields  (5.8±  1 .8)x  10” cm' 

molecule"'  s"'.  The  joint  uncertainty  limits  thus  do  not  quite 


overlap.  The  value  of  a?  =0.5  in  Eq.  (3)  was  chosen  for  uni¬ 
formity  with  the  work  on  the  .v'  and  s~  metals, '  “  and  may  not 
be  the  best  choice  here.  Also  extrapolation  over  a  large  range 
of  7  ‘ .  i.e.,  from  ( 1 .9  to  3.4) x  10  ’  is  a  .somewhat  uncertain 
procedure.  1  here  is  no  reason  to  doubt  either  data  set.  In¬ 
deed,  we  recently"'’  made  measurements  on  the  Cr  +  N^O  re¬ 
action  in  a  MHTP  (metals  high-temperature  photochemistry 
reactor)  from  278  to  1150  K,  which  showed  excellent  agree¬ 
ment  with  the  static  cell  results  at  298  and  348  K,'^  as  well 
as  with  a  shock  tube  study‘s  from  1130  to  2570  K.  In  the 
same  work  we  reported  similar  agreement  between  a  MHTP 
study  of  the  Cr+HCl  reaction  from  811  to  1414  K  and  HT- 
FFR  measurements  from  880  to  1449  K.  As  the  various  tech¬ 
niques  for  studies  of  elementary  metal  atom  reactions  are 
showing  good  agreement,  .see  also  Refs.  5  and  18,  it  appears 
justified  to  combine  the  two  sets  of  results  for  the  present 
reaction.  Using  the  fitting  procedures  given  above  leads  to 

A'(296-1030  K) 

=  5.8X  10"'^7/K)'-''^ 

Xexp( -0.752//?7)cm^  molecule"'  s"'  (4) 

with  (7^  =  3. 19  A~,  cr^„  =  “0.428  A,  cr^£;  =  2.86  AE, 
(;'  =  7.86X  I0'\  -0.382  £,  and  oi=2.71  are 

yielding  7:20*  precision  limits  from  ±3%  to  ±21%  and  ac¬ 
curacy  limits  from  ±23%  to  ±31%.  Figure  2  shows  that  this 
expression  leads,  in  the  temperature  range  of  the  present 
measurements,  to  a  fit  nearly  as  good  as  that  resulting  from 
Eq.  (3). 


III.  THEORY 

A.  Activation  barriers  of  Z(s*p’)+N20(X 
reactions  using  semiempirical  configuration 
interaction  (SECI)  theory 

Following  the  treatment  of  Ref.  1,  the  wave  function  i/» 
of  the  activated  complex  (NNOZ)*  can  be  taken  as  the  linear 
superposition  of  three  wave  functions  (Ai  •  corre- 

spending  to  the  properties  of  three  resonating  structures. 
The.se  result,  respectively,  from  the  covalent  ground-ground 
state  interaction,  a  covalent  interaction  where  a  valence  s 
electron  is  promoted  to  a  state,  and  an  ionic  interaction 
between  the  positively  charged  atom  and  negatively  charged 
NoO.  Thus, 

<//=  c  ]  +  <-*2  ^  1 

where  (A  =  1,2,3)  are  the  expansion  coefficients.  The 
height  of  the  barrier  is 

According  to  the  variation  principle, the  wave  function  4/ 
for  the  system  is  that  for  which  E  is  at  a  minimum.  This 
condition  can  be  expres.sed  as  a  set  of  simultaneous  linear 
equations, 

2  Q(/y;,-A/;£)  =  0,  /=  1,2,3.  (7) 

k  1 
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T.\Bl,h  I.  .Suinniary  nf  rate  coetMcienJ  measurements  ot  Al  f  N-.O.'^ 


Reactidn 

F 

Initial  Al 

/one  length 

P 

l.v/l 

(arbitrary 

absorption 

V 

7 

(cmi 

1  mhur) 

I  lO' '  cm  ') 

( 1  {)'“  cm  ') 

units) 

P/r) 

(ms  ' ) 

(K) 

20 

18.2 

2.1 

3.7 

84 

28 

623 

20 

18.2 

2.1 

3.6 

69 

28 
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20 

18.2 

2.0 

3.5 

62 

29 

646 

20 
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3.7 

3.8 

52 

27 

672 

20 

.M.6 

3.7 

3.9 

51 

26 

670 

10 

.^4.5 

3.9 

4.0 

70 

26 

646 

10 

16.6 

1.9 

5.0 

17 

21 

629 

20 

44.9 

4.4 

2.6 

47 

17 

745 

20 

44.9  . 

4.2 

2.6 

57 

18 

765 

20 

44.7 

4.2 

2.5 

42 

18 

771 

20 

.M.i 

3.9 

2.5 

46 

18 

640 

20 

M.2 

3.8 

4.1 

72 

19 

655 

20 

20.3 

2.2 

2.9 

102 

35 

665 

10 

47.9 

4.5 

6.0 

49 

17 

767 

10 

47.9 

4.5 

6.0 

52 

17 

771 

10 

28.3 

2-6 

4.4 

68 

29 

776 

!() 

54.4 

5.0 

6.0 

51 

17 

796 

10 

52.9 

4.4 

4.1 

68 

18 

864 

20 

48.7 

3.8 

4.6 

47 

22 

928 

20 

48.4 

3.8 

3.2 

29 

22 

927 

20 

45.7 

3.6 

3.4 

49 

30 

921 

20 

49.3 

3.9 

4.2 

38 

30 

915 

20 

19.4 

1.5 

2.6 

51 

40 

938 

20 

23.7 

1.8 

3.8 

57 

47 

931 

20 

22.7 

1.8 

3.0 

40 

60 

921 

20 

22.6 

2.1 

3.1 

34 

33 

784 

20 

22.6 

2.1 

3.9 

36 

33 

789 

20 

25.8 

2.4 

4.5 

31 

34 

792 

20 

31.0 

2.8 

4.2 

26 

36 

796 

10 

20.0 

1.9 

4.9 

43 

31 

774 

10 

20.0 

1.9 

5.7 

34 

31 

775 

10 

34.2 

3.2 

7.1 

21 

25 

775 

20 

21.7 

3.0 

4.7 

37 

27 

516 

20 

28.9 

3.5 

5.8 

44 

22 

60! 

20 

17.3 

1.9 

12.9 

57 

47 

671 

20 

(7.3 

1.9 

12.9 

56 

47 

671 

10 

17.3 

1.9 

13.2 

62 

46 

656 

10 

16.0 

1.8 

13.6 

68 

45 

656 

10 

20.8 

2.2 

12.2 

63 

50 

675 

20 

16.1 

1.4 

9.1 

24 

57 

826 

20 

16.1 

1.4 

10.8 

48 

57 

826 

20 

14.7 

1.3 

11.1 

46 

55 

834 

20 

14.7 

1.3 

10.9 

32 

56 

847 

10 

14.7 

1.3 

11.4 

34 

54 

813 

10 

14.7 

1.3 

13.1 

39 

54 

812 

10 

14.7 

1.2 

12.0 

56 

58 

881 

10 

14.7 

1.2 

11.8 

47 

60 

896 

10 

14.8 

1.2 

il.8 

44 

60 

906 

10 

14.8 

1.2 

11.6 

28 

61 

917 

20 

16.1 

1.3 

9.7 

52 

63 

913 

20 

16.1 

1.3 

11.3 

69 

63 

20 

19.8 

1.6 

9.3 

35 

66 

901 

20 

19.9 

1.6 

9.3 

26 

66 

900 

20 

26.0 

2.5 

4.6 

9 

26 

743 

20 

39.5 

3.8 

4.8 

25 

21 

751 

20 

39.5 

3.8 

4.8 

22 

21 

752 

20 

16.4 

1.7 

8.2 

21 

50 

681 

20 

16.4 

1.7 

10.9 

25 

51 

683 

10 

16.3 

1.9 

17.5 

24 

47 

632 

10 

17.2 

2.2 

13.9 

16 

49 

568 

10 

17.3 

2.2 

27.7 

31 

49 

574 

20 

17.2 

2.1 

27.3 

28 

50 

580 

20 

17.2 

2.1 

16.4 

21 

50 

582 

20 

16.8 

1.3 

8.1 

15 

51 

963 

20 

20.3 

1.5 

7.4 

14 

56 

991 

20 

20.3 

1.5 

7.4 

15 

56 

992 

20 

19.2 

1.4 

7.1 

26 

77 

996 

10 

19.5 

1.5 

11.1 

15 

74 

968 

10 

19.5 

(.5 

18.4 

20 

74 

967 

20 

32.8 

3.8 

8.4 

16 

19 

627 

20 

37.9 

2.7 

4.6 

15 

26 

1024 

20 

37.9 

2.7 

4.6 

17 

26 

1029 

The  measurements  are  reported  in  the  sequence  in  which  they  were  obtained. 
^Should  be  read  as  (4.49±0.42)X  10’“, 


k  z 

Icin' molecule  '  s  ') 


4.49i0.42  (-11)^ 
4.33±0.43  (-11) 
4.33 ±0.40  (-11) 
4.73±0.47  (-11) 
4,35±0,40  (-11) 
4.25±0.46  (-11) 
4.75 ±0.65  (-11) 
5.1I±0.40(-fl) 
5.41  ±0.45  (-1!) 
5.91  ±0.46  (-11) 
4.39±0.45  (-11) 
4.I2±0.36  (-11) 
3.46±0.37  (-11) 
5.29±0.39  (-11) 
5.44±0.45  (-11) 
5.59 ±0.50  (-11) 
6.3 1  ±0.54  (-11) 
7.52±0.75  (-1!) 
6.76±0.53  (-11) 
7.I8±0.48  (-11) 
6.86±0.62  (-!!) 
7.28±0.54  (-11) 
6.65 ±0.75  (-11) 
6.71  ±0.80  (-11) 
6.00±0.62  (-11) 
5.77±0.57  (-11) 
4.22±0.35  (-11) 
5.07 ±0.66  (-11) 
5.01  ±0.51  (-11) 
4.37±0.55  (-11) 
3.72 ±0.44  (-11) 
3.53 ±0.45  (-11) 
2.40±0.26  (-11) 
2.52±0.25  (-11) 
3.19±0.38  (-11) 
2.93±0.32  (-11) 
3.58±0.41  (-11) 
3.61  ±0.59  (-11) 
4.96±0.47  (-11) 
4.57±0.49  (-11) 
4.56±0.48  (-11) 
6.05 ±0.66  (-11) 
5.97±0.76  (-11) 
4.17±0.57  (-11) 
4.53 ±0.74  (-11) 
5.95±0.87  (-11) 
6.53 ±0.74  (-11) 
4.94±0.66  (-11) 
6.I6±0.78  (-11) 
7.58±0.82  (-11) 
7.22±0.75  (-11) 
7.55 ±0.68  (-11) 
7.94±0.70  (-11) 
3.07±0.24  (-II) 
4.09±0.40  (-11) 
3.77±0.28  {-ID 
3.91±0.41  (-11) 

3.5 1  ±0.35  (-11) 
3.35 ±0.37  (-11) 
3.01  ±0.34  (-11) 
2.50±().30  (-11) 

2.3 1  ±0.24  (-11) 
2.66±0.31  (-11) 
5.63±().58  (-11) 
6.36±0.55  (-11) 
6.63 ±0.59  (-11) 
7.71  ±0.70  (-11) 
6,()3±0.65  (-11) 
5.57±().55  (-11) 
3.46±0.29  (-11) 
5.62±0.39  (-11) 
5.75 ±0.47  (-11) 
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T(K) 


FIG.  1.  Arrhenius  pk>t  nt  the  rale  La>etticient  data  Al^N.O  reaction. 
;  j _ LIP  used  as  the  meniloring  technique.  .ARAS  used  as  the  monitor¬ 

ing  technique. 


T(K) 


FIG.  Z.  Arrhenius  plot  of  all  the  rate  coefhcienl  data  tor  Al-^-N.O  reaction: 
_ Present  work.  □ — Parnis,  Mitchell.  Kanigan,  and  Hackett  (Ret.  12). 


Here  H represents  the  interaction-energy  integrals  and  A//, 
the  overlap  integrals,  which  are  defined  as 

W,*=|  ar.  li,=  j  ctr.  (8) 

The  overlap  integrals  are  approximately  equal  to  zero.  A/^=0 
for  because  each  of  the  three  resonating  structures  are 

distinctly  different  from  the  others.  This,  along  with  normal¬ 
ization  of  Eq.  (5).  results  in  Aa  =  /i/'f  'A* 


I 

where  the  cl  represent  the  fraction  of  character  contributed 
by  the  respective  structures  to  the  activated  comple.x.  Simi- 
larly,  Eq.  (7)  reduces  to  the  secular  equation 


W..  =  PE-Hi/(iV')/2.  ("> 

The  third  structure  results  from  the  ionic  interaction  of 
Z'"(.s-)  +  N.O‘{-n).  which  involves  the  attractive  Coulomb 
and  the  repulsive  Born  potential  between  these  species.  The 
energy  of  the  barrier  resulting  from  this  structure  is  evaluated 
as  before" 


where  IP  is  the  ionization  potential  of  Z.  E,^  is  the  vertical 
electron  affinity  of  N.O,  e  is  the  elementary  charge.  /?„  is  the 
sum  of  the  univalent  radii  of  Z’^  and  O  .  and  S  is  the  average 
Born  e.xponent  of  Z "  and  O'.  Values  for  these  properties  and 
PE  along  with  their  sources  are  given  in  Table  II. 

The  interaction-energy  integrals  are  estimated  to  be 
equal  to  those  for  the  alkali,  alkaline  earth,  and  transition 


,i-£ 

W,: 

Wl2 

W,, 

H.y 


=  0. 


(10) 


w,,  is  the  energy  of  the  Z(s-p')  +  n,OiX  'S^)  interac¬ 
tion.  The  barrier  height  for  this  ^tructure  qis'p')  equals 
W|  I .  For  each  series  of  homologous  reactions  q  is  constant; 
it  is  the  one  quantity  to  be  determined  from  experiment  in 
these  treatments.  The  second  ^lructure  results  from  the  at¬ 
tractive  covalent  interaction  of  Zi.v'p'i  with  NiO  where  spin 
conservation  requires  N.O  to  be  in  a  triplet  state,  whereas  tor 
the  ,9'  and  .v"  metals  N.O  is  in  a  kinglet  state  in  the  second 
structure.  The  barrier  to  form  this  structure  includes  the  pro¬ 
motion  energy  PE  of  Z(s-p')  to  the  Zlx'/)")  state  and  the 
e.xcitation  energy  of  N.OiX  A ' )  to  the  hrst  N.O(  I' )  state, 
minus  the  amount  by  which  the  interaction  is  attractive.  Fol¬ 
lowing  the  previous  work'  "  this  term  is  estimated  to  be 
q(s-p')l2.  Thus, 


TABLE  11.  Physical  properties 
13  atom  reactions  with  NnO. 

used  to  cal 

Iculale  energy  integrals 

for  group 

PE 

IP-EA 

Rp 

Z 

(kJ  mol  '  i 

ikJ  mol  ‘ 

)  (  to  "  cm) 

B 

U5'' 

10!  S’’ 

2.2.T 

h,(F 

A1 

347 

705 

2.53 

7.0 

Ga 

45. S 

706 

2.5H 

S.O 

In 

410 

776 

S,5 

T1 

54! 

S07 

3.05 

0  5 

The  values  correspond  to  the  lowest  v-p  promotion  enerev  obtained  trom 
the  atomic  energy  levels  in  Ret.  24. 

"Values  are  tihlamed  fnim  the  lonizatum  potentials  m  Ret.  >>  and  the  ver- 
„cal  electron  .ilhnilv  ol  N-O.  -217  U  mol  which  corresponds  to  the 
N4.)  negative-ion-state  "11  to'in  Ret.  26.  n  ■  -.-i  tk 

Univalent  radii  are  derived  from  the  etfcLTive  ionic  radii  in  Ret.  .e 

univalent  radius  of  O  is  taken  to  be  1.76X10  vm.  .  „ 

^'The  average  Born  exponent  equals  ((5, +7 i/2  where  the  value  ot  /  is  t  e 
Born  exponent  of  O  and  the  (5,  value  is  that  ot  Z'  taken  trom  Ret. 
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TABLE  III.  Compan.son  of  e.xp^iri mental  and  calculated  activation  barriers  tor  group  13  atom  reactions  with 


N.O. 


Z 

<-■1 

(fraction) 

(fraction) 

(fraction) 

4 

(cm^  molecule"’  s" ‘ 

)  (kJ  mol  ' 

F 

)  (kJ  mol " ’) 

(K) 

B 

0.85 

0.09 

0.06 

8.4X10"’^ 

21.3 

22.1'' 

300 

Ref.  20'' 

A1 

0.79 

0.09 

0.12 

4.6X10"'^“ 

6.5 

6.5 

520-1030 

Thi.s  work 

Ga 

0.81 

0.06 

0.13 

5.9X10“'- 

14.3 

14.4" 

296 

Ref.  21'' 

In 

0.81 

0.07 

0.12 

6.3X10“'^ 

13.4 

T1 

0.87 

0.04 

0.09 

7.1X10“'^ 

26.7 

‘‘Experimental  value  of  the  AI  +  N2O  reaction;  A  coefficients  for  all  other  reactions  are  calculated  using  simple 
collision  theory. 

^is  £:„p,  value  is  obtained  from  the  300  K  rate  coefficient  and  the  A  coefficient  of  column  4. 

"This  £„p,  value  is  obtained  from  the  296  K  rate  coefficient  and  the  A  coefficient  of  column  4. 

‘^Reference  from  which  the  experimental  data  were  obtained. 


metal  atom  reactions.”  Thus  H ,2  and  //,3  are  equal  to  /3  and 
H^^  is  equal  to  /3/2,  where  ^=-107.9  kJ  mol 

6.  Results  and  discussion 

The  value  of  q{s"p^)  is  chosen  such  that  the  height  of 
the  barrier  E  of  Eq.  (10)  equals  the  experimental  activation 
barrier  E^^pt  of  the  AI  +  N2O  reaction  (Sec.  II  B).  From  Eq. 
(3)  this  yields  <y(5  V)=85.7  kJ  mor'.  Using  this  value,  and 
Eqs.  (7)  and  (9),  the  activation  barriers  E  and  the  fraction  of 
character  contributed  by  each  structure  can  now  be  calcu^ 
lated  for  the  other  boron  group  atoms.  The  results  of  these 
calculations  are  given  in  Table  III.  Table  III  also  shows  the 
comparison  of  the  calculated  activation  barriers  to  the  de¬ 
rived  experimental  activation  barriers  E^,p^.  For  B  and  Ga, 
these  experimental  activation  barriers  are  determined  by  set¬ 
ting  the  rate  coefficients  found  in  room  temperature 
studies”'^"'  equal  to  exp(-- E/EE)  in  which  the  A  are 


FIG.  3.  Plot  of  the  experimental  and  calculated  activation  barriers  (crosses) 
of  Z^NLO— ZO  +  N2  reactions  vs  the  sums  of  the  s-p  promotion  energies 
and  ionization  potentials  of  Z.  .Mkali  metal  (open  triangles),  alkaline  earth, 
and  transition  metal  (open  circles),  and  boron  group  (open  squares)  atom 
reactions  are  distinguished. 


taken  as  the  values  given  by  simple  collision  theory”^ 
Rl(^7rkg/pL)^'^.  Here,  kg  is  the  Boltzmann  constant  and  /x 
the  reduced  mass.  The  predicted  barriers  for  B  and  Ga  are 
within  0.7  kJ  mol”’  of  the  derived  experimental  values,  us¬ 
ing  the  lowest  s-p  promotion  energies.  Using  instead  the 
activation  barrier  from  the  combined  experimental  rate  coef¬ 
ficient  expression,  Eq.  (4),  yields  values  within  1.2  kJ  mol"' 
of  the  predicted  values.  This  confirms  that  varying  the  value 
of  n  has  little  influence  on  the  accuracy  of  the  predictions. 
The  Eg^pt  for  In  and  TI  +  N2O  could  not  be  estimated  because 
no  experimental  data  are  available  for  these  reactions. 

In  the  present  treatment,  Sec.  Ill  A,  we  assumed  the 
lowest  s-p  promotion.  The  activation  barriers  were  also 
computed  using  the  lowest  p-s  promotion  energies  and  the 
average  of  the  promotion  energies  \_{s-p)  +  {p-s)]l2.  In  all 
cases  the  calculated  activation  barriers  were  too  high  for  bo¬ 
ron  and  too  low  for  gallium,  indicating  that  the  lowest  p-s 
promotion  energies  do  not  correlate  with  the  activation  bar¬ 
rier  of  these  atoms. 

In  Fig.  3  the  experimental  and  calculated  activation  bar¬ 
riers  of  the  present  work  are  compared  to  those  for  the  5  ‘  and 
metals  atom  N2O  reactions.''^  As  can  be  seen,  good  agree¬ 
ment  between  theory  and  experiment  is  obtained  for  ail  three 
groups.”^  Each  group  shows  a  linear  increase  of  the  activa¬ 
tion  barrier  with  the  sums  of  the  ionization  potentials  and 
s-p  promotion  energies,  which  properties  determine  the  en¬ 
ergies  of  the  resonating  structures,  cf.  Eqs.  (11)  and  (12).  As 
the  magnitude  of  IP  +  PE  increases,  the  resonance  interaction 
between  these  structures  and  the  activated  complex  weakens, 
thereby  increasing  the  barrier  E. 

In  summary,  the  semiempirical  treatment,  applied  to  5' 
and  s"  metal  atom  reactions  with  N2O,  has  been  extended  to 
atoms  by  modifications  to  the  treatment  and  use  of  the 
experimental  barrier  from  one  5”p'  atom  reaction,  namely 
AI-fN^O.  Comparison  to  the  experimental  data  available  in 
the  literature  indicates  good  agreement  between  theory  and 
experiment. 
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The  AlO  +  O2  Reaction  System  Over  A  Wide  Temperature  Range. 
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High-Temperature  Reaction  Kinetics  Laboratory,  The  Isermann  Department  of  Chemical 
Engineering,  Rensselaer  Polytechnic  Institute,  Troy,  NY  12180-3590 
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The  kinetics  of  the  AlO  +  O2  "action  system  has  been  studied  in  a  high-temperature  fast-flow 
reactor  (HTFFR)  from  7  to  100  mbar.  The  relative  AlO  concentrations  were  monitored  by  laser- 
induced  fluorescence  at  the  and  C2n-X2l  transitions.  Below  about  1000  K  the  AlO 

consumption  is  pressure-dependent,  i.e.  leading  initially  to  AIO3,  with  log  k(305-1010  K)  = 

-25.36  -  1.69  log  (T  /  K)  cm^  molecule'^  s'^  as  obtained  by  weighted  linear  regression.  At  higher 
temperatures  the  system  becomes  pressure-independent  with  AIO2  and  O  as  the  likely  products.  For 
this  process  k(  1200-1690  K)  =  7.7  x  lO'lO  exp  (-10008  K  /  T)  cm3  molecule-1  5-1.  The  accuracy 
estimates  are  discussed  in  the  text.  The  reactions  are  suggested  to  proceed  through  an  OAIOO  addition 
complex  of  similar  structure  as  identified  in  solid  matrices. 


INTRODUCTION 

Reactions  of  aluminum  species  play  a  role  in  many  high-temperature  environments.  We 
have  studied  a  large  number  of  Al,  AlCl  and  AlO  reactions  over  wide  temperature  ranges,  by  use  of 
the  H'I'FFR  (high-temperature  fast-flow  reactor)  technique.^-*  Two  early  studies  of  the 
AlO  +  O2  system,  at  respectively  1400  and  300  K,9d0  suggested  a  temperature-independent  rate 
coefficient,  attributed  to  abstraction 

AlO  +  O2  — ^  AIO2  +0  (1) 

This  interpretation,  which  indicates  no  activation  energy,  is  in  conflict  with  some  subsequent 
information,  such  as  the  currently  recommended  endothermicityl  ^  for  reaction  (1),  AH°  =  96  +  40 
kJ  mol"'.  Recently,  we  showed  that  several  reactions,  which  were  considered  in  the  literature  to 
proceed  by  simple  abstraction,  i.e.  AlO  +  CH4  and  Cr  +  O2,  to  be  more  complex .6-12  The  earlier 
AlO  +  O2  measurements  could  have  similarly  hidden  an  addition  process,  as  they  showed 
considerable  scatter.  Since  then  the  HTFFR  technique  has  been  much  improved.5d3  a  re¬ 
investigation  of  the  AIO/O2  system  thus  was  indicated.  We  demonstrate  here  that  there  are  in  fact 
two  reactions,  reaction  (1)  which  under  the  present  conditions  dominates  above  1200  K,  and  an 
addition  reaction 


AlO  +  O2  +  Ar  — ^  AIO3  +  At 
which  dominates  below  1000  K. 


(2) 


TECHNIQUE 

The  HTFFR  facility  used  in  this  work  and  the  data  handling  procedures  have  been  previously 
described.^d4  Briefly,  a  vertical  mullite  (McDanel  MV-30)  reaction  tube  (60  cm  long,  2.2  cm  i.d.), 
was  heated  to  the  desired  temperature  by  resistively  heated  SiC  rods  inside  an  insulated,  water-cooled 


vacuum  housing.  AlO  was  produced  by  several  methods.  At  room  temperature  and  for  a  few  checks  at 
547  K,  production  was  achieved  by  passing  Ar  with  approximately  100  ppm  of  trimethyl  aluminum 
(TMA)  and  1  ppm  of  O2  through  a  microwave  discharge.  This  technique  was  limited  to  reaction  zone 
pressures  from  7  to  28  mbar  and  could  not  be  used  at  higher  temperatures.  For  the  other  experiments 
A1  vapor  was  produced  by  resistively  heating  an  Al-wetted  W  coil  in  a  flow  of  Ar  carrier  gas.  A  trace  of 
oxidizer,  0.5%  N2O  or  5%  O2  in  Ar,  at  c  4  x  10-3%  of  the  main  Ar  flow,  was  passed  through  a  Pt 

tube  to  a  location  just  downstream  from  the  vapor  source  to  rapidly  convert  A1  to  AlO.^’^’^^’^^  This 
production  technique  limited  observations  at  547  K  to  P  <  32  mbar  due  to  the  condensation  of  A1  and 
AlO  on  the  reactor  walls,  but  at  higher  ternperamres  allowed  observations  up  to  100  mbar,  i.e.  over  the 
full  normal  HTFFR  pressure  range.  For  the  T  >  850  K  experiments  only  the  5%  O2  was  used. 

Downstream  from  the  AlO  production  zone,  02/Ar  mixtures  at  flow  rates  between  1%  and  8% 
of  the  main  Ar  flow  rate  were  introduced  through  a  movable  inlet  at  10  or  20  cm  upstream  from  the 
observation  windows.  There,  relative  AlO  concentrations  were  measured  by  laser-induced  fluorescence 
(LEF).  A  pulsed  Lambda  Physik  EMG  101  excimer  /  FL  2002  dye  laser  in  combination  with  a  KDP 
doubling  crystal  was  employed.  Just  before  the  reactor  window,  the  laser  beam  passed  through  two 
light  baffles, >7  made  of  metallic  concentric  brass  cones  with  4  and  2.5  mm  diam.  exits  respectively,  to 
reduce  scattered  light  in  the  reactor.  Two  transitions,  B2E-X2I  and  were  used.  For  the 

former,  AlO  is  pumped  on  the  464.8  nm  (1,0)  band  and  the  fluorescence  is  observed  through  a  482nm 
(20  nm  full  width  at  half  maximum)  interference  filter,  i.e.  mainly  at  the  486.6  nm  ( 1,1)  band. The 
B3S-X2Z  system  could  only  be  used  up  to  1200  K  because  of  interference  by  background  radiation 
from  the  reactor  walls.  For  measurements  above  1200  K  and  for  a  comparison  near  826  K  the  weaker 
LIF  from  the  C^n-X^I  (0.0)  transition  at  302.2  nm'^  was  used  in  combination  with  a  301nm  ( 1 1  nm 
full  width  at  half  maximum)  interference  filter.  The  fluorescence  was  detected  by  an  EMI  98 13QA 
photomultiplier  tube,  connected  to  a  Data  Precision  Analogic  6000/620  100  MHz  transient  digitizer. 
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The  gases  used  were  Ar  (99.998%)  from  the  liquid  (Linde),  O,  (99.6%)  from  Linde,  O, 
(99.98%)  from  Matheson  and  5%  (99.99%)  in  Ar  (99.995%)  from  Scott.  The  rate  coefficients 

obtained  with  these  O2  sources  were  indistinguishable.  The  5%  O2  and  0.5%  O2  (99.99%)  m  Ar 
(99.995%)  from  Matheson,  0.5%  N2O  (99.99%)  in  Ar  (99.995%)  from  Matheson,  and  TMA  (98%) 
from  the  liquid  TMA  (Johnson  Matthey)  were  used  for  AlO  production. 


Measurements  were  made  under  pseudo-first-order  conditions,  [AlO]  « [O2],  in  the  stauonary 
inlet  mode.^^  Rate  coefficients  ki,  in  second  order  units,  for  a  given  temperature,  pressure,  and  average 
velocity  were  obtained  by  using  five  O2  concentrations  providing  variation  by  a  factor  of  about  5.  ki 
and  CTki  were  calculated  by  applying  a  weighted  linear  regression^O  to  plots  of  In  [AlOlreiative  ^s.  [O2]. 
These  yielded  straight  lines  with  slopes  equal  to  -kit,  where  t  is  reaction  time. 


RESULTS 

The  measurements  are  divided  into  two  sets.  In  one,  the  measurements  were  made  at  varying 
pressures  at  temperatures  near  305,  547, 826.  or  1010  K  to  evaluate  the  temperature  dependence  of  the 
ptessute-dependent  reaction.  In  the  other,  T  >  1200  K,  where  pressure-independent  rate  coefficients 
were  obtained,  the  temperature  was  randomly  varied  up  to  1690  K. 

305-1010  K  MEASUREMENTS 

The  ki  measurements  obtained  in  this  range,  and  the  conditions  under  which  they  were  obtained, 
are  assembled  in  Table  I  and  show  a  pressure  dependence.  At  each  temperature,  to  determine  if  the 
data  were  dependent  on  parameters  other  than  [M],  plots  of  the  residuals  [k;  -  k(M)]/k(M)  versus  an 
individual  parameter  were  examined.  Here,  k(M)  is  the  fitting  expression  for  that  particular  temperature, 
eqs.  3,  4,  or  5  below.  These  plots  showed  that  the  third-order  rate  coefficients  are  independent  of  initial 
fluorescence  intensity  F  (a  measure  of  initial  [AlO] ),  average  gas  velocity  T,  and  observed  reaction  zone 


length. 
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There  were  slight  variations  in  temperature  for  each  group  of  measurements.  The  standard 
deviation  of  these  were  computed2i  for  each  data  set  and  combined,  by  a  propagation  of  error 
method, 22  with  an  estimated  maximum  tTT  =  ±  2%  T  for  systematic  errors.  The  resulting  standard 

deviations  for  each  average  temperature  are  shown  in  Table  II. 


The  305  K  measurements  are  near  the  high-pressure  limit  as  evidenced  by  the  leveling  off  of  the 
kj  vs.  [M]  plot.  Fig.  1.  These  data  were  fitted  using  Marquardf  s  method,23  with  a  fall-off  expression 

fromTroe;24 


logk(M)  =  log! 


k2.o[M] 


+  k2,o[M]/k2 


log  Fc 


l-i-[log(k2,o[M]/k2,oo)l 


(3) 


Here  k2,o  is  the  third-order  rate  coefficient  for  reaction  (2),  kj...  is  the  high-pressure  second-order  rate 
coefficient,  and  Fc  is  a  broadening  factor.  Eq.  (3)  y.elds  kj.o  =  5. 1  x  10-30  c„6  ntolecule'^,  s' ' . 
k,.=  l.l  X  10-12  cm3  molecule-1  s-',  and  Fc  =  0.752.  To  estimate  Okj.o  a  "<>”'>"0“  f"”  “‘  *0 
was  made  with  only  the  first  term  of  eq.  3.  From  this  the  elements  of  the  covariance  matrix  were 
calculated25  and<Tk2.o  was  found  to  be  0.34  k2,0r  Table  II. 

At  547  and  826  K  linear  k,  vs.  [M]  plots  were  obtained.  Figs  2  and  3.  The  kj  were  fitted  with 
an  expression  purely  proportional  to  pressure, 

k(M)=  k2.o  [M] 

From  the  fits23  the  covariance  matrices  were  calculated23  and  the  Oki.o  obtained.  Table  II.  The  .s47  K 
measurements  allowed  comparison  of  the  various  AlO  production  techniques.  As  shown  in  Fig.  2  no 
infiuence  of  the  method  used  is  evident.  At  826  K  higher  pressures  could  be  studied  than  at  .547  K.  due 


to  the  improved  efficiency  of  the  production  technique.  The  larger  [AlO]  at  this  temperature  also  made  it 
convenient  to  compare  the  LIF  results  from  the  transition  to  those  of  the  C^n-X^I 

transition.  This  is  done  in  Fig.  3,  which  shows  the  independence  of  the  transition  used. 

The  data  at  1010  K  display  linear  behavior  with  an  intercept  indicative  of  a  pressure-independent 
channel,  Fig.  4.  These  data  were  fitted  with  a  simple  linear  expression 

k(M)  =  ki  +  k2,o  [M]  (5) 

using  a  linear  least-squares  algorithm.^®  In  this  eq.  ki  is  the  second-order  rate  coefficient.  This 
algorithm  also  yielded  Ok2,o 

The  temperature  dependence  of  the  third-order  rate  coefficients  is  shown  in  Fig.  5.  The  rate 
coefficients  decrease  with  increasing  temperature,  typical  behavior  for  an  association  reaction.  The  data 
were  fitted  with  a  log  k2,o  =  log  A2,o  +  n  log(T  /K)  expression  using  a  linear  least  squares  fit^O  with  the 
transformation  Cj  =  2.303  Toiog  j  which  yields 

log  k2.o  (305  -  1010  K)  =  -25.36  -  1.69  log  (T  /  K)  cm^  molecule-2  g-l.  (6) 

The  variance  and  covariance  matrix  elements25  are  (T2iogA2.0  =  1-45,  cTiogA2.0  n  =  -0.508,  and 
ct2„  =  0. 178.  Using  these  variances  and  covariance  the  2ok2,o  Precision  limits  are  calculated  to  vary 
from  ±76%  at  305  K,  decreasing  to  a  minimum  of  ±22%  at  700  K,  and  increasing  again  to  ±37%  at 
1010  K.  Using  a  ±10%  uncertainty  in  the  reactor  flow  profile,^^'"^  along  with  a  conservative  ±20% 
for  possible  systematic  errors,  yields  confidence  intervals  varying  from  ±79%  to  ±31%. 
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1200  -  1690  K  MEASUREMENTS 

The  1200  K  to  1690  K  data  are  assembled  in  Table  HI.  Res, dual  plots  showed  no  dependence  on 
pressum,  initial  fluorescence  intensity  F,  average  gas  velocity  7,  and  reaction  zone  length.  As  a 
precautton.  in  deterrmning  the  parameters  of  the  lc,(T)  e>tpression  below,  only  measurements  a,  low 
enough  pressures  to  insure  less  than  20%  contribution  from  dre  extrapolated  third-order  reaction  eq.  6 
were  used.  These  measurements  are  shown  in  Fig.  6  and  show  a  strong  positive  temperature  dependence 
with  no  significant  curvature.  Applying  a  linear  regression  fit^S  to 
ki(T)  =  Aexp(  -E  K/  T)  yields 


k, (1200-1690  K)  =  7.7x  10-10  exp(-10008K/T)  cm’ molecule' >  s 


(7) 


A  nearly  identical  expression  was  obtained  by  (i)  fltting  all  the  data  and  (ii)  by  subtractmg  the 
termolecular  component  from  all  the  data  and  fining  the  results.  In  both  cases  dre  acdvation  energy  w^ 
withrn  6%  of  the  value  shown  in  eq.  7.  Equation  7  has  variances  and  covariance’’  of  tr,’  =  0. 182  A’, 
(T.,£  =  2.54  X  10’  A,  and  =  3.59  x  10’,  which  yields  precision  limits  varying  from  +9%  to  ±17% 
depending  upon  temperature.  The  resulting  estimated  accuracy  limit  is  ±26%  when  the  same  systematic 
errors  as  above  are  considered.  The  k,  from  the  1010  K  measurements  (above)  was  no,  used  as  the 
associated  uncertainly  was  considered  to  be  too  large. 


DISCUSSION 

The  results  of  die  negative  temperaniie  dependence  of  reaction  (2)  and  the  positive  dependence 
of  reaction  ( 1)  is  that  the  AlO  consumption  rate  coefficients,  calculated  on  a  second  order  basis,  pass 
through  a  minimum  which  occurs  al  around  800  -  1000  K  at  the  pressures  used.  Comparison  of  the  k, 
data  at  300  K  (Table  I)  to  those  at  1400  K  (Table  ID)  shows  little  difference.  Thus  the  temperature 
independence  indicated  for  reaclion  ( 1 )  in  the  previous  work'd  „as  based  on  a  coincidence.  The  actual 
measurements  from  the  present  and  previous’.  ">  work  are  in  reasonable  agreement  as  can  be  seen  in 
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Fiss.  5  and  6.  For  the  comparison  in  Fig.  5  the  individual  second-order  rate  coefficients  from  ref.  10 
were  divided  by  the  pressure  and  averaged.  There  now  is  also  no  conflict  between  the  endothermicity 
for  reaction  (1)  of  96  ±  40  kJ  mol"*  and  the  activation  energy  of  83  ±  10  kJ  moP^  obtained  from  eq.  7. 

The  type  of  results  obtained  here  can  be  fitted  to  a  general  scheme  of  the  kind: 

AlO  +  O2  AIO3  AIO3 
!■ - ■»  AIO2  +  O 

,  where  the  dotted  lines  indicate  the  less  likely  channels.  The  path  leading  to  AIO3  appears  a  straight 
forward  termolecular  mechanism.  The  chemically  activated  intermediate  AIO3  is  also  likely  to  be  the 

precursor  of  the  AIO2  +  O  reaction  products  of  the  high-temperature  second-order  process,  as  the  pre- 
exponential  of  the  ki(T)  expression,  7.7  x  10'^9  molecule"^  s'l,  appears  rather  too  large  for  a 
diatomic  radical  -  diatomic  radical  direct  abstraction  reaction.27  Certainty  with  regard  to  the  mechanism 
of  the  relative  importance  of  direct  and  indirect  bimolecular  processes  is,  however,  generally  difficult  to 
establish.-^ 

The  large  pre-exponential  of  ki(T)  also  suggests  that  the  AIO3  complex  once  formed  will  not 
significantly  re-dissociate  to  the  reactants  AlO  and  02-  This  in  contrast  to  the  exothermic  iso-electronic 
BO  -I-  O2  reaction  system  for  which  no  pressure  dependence  and  a  negative  activation  energy  was 
found. This  suggests^^’^9  ^  reaction  invob-  ng  an  activated  complex  which  is  not  significantly 
pressure  stabilized.  The  mechanism  for  that  reaction  can  be  summarized  by 

BO  +  O2  ^  BO3  — >  BO2  +  O 


There  have  apparently  not  been  any  direct  gas-phase  observations  of  AIO3  species.  However, 

the  A1  +  O2  and  A1  +  O3  reactions  when  carried  out  in  solid  matrices  have  been  shown  to  yield  a 

product3h32  OAIO2  that  is  also  the  likely  product  for  AlO,  O2  addition.  The  matrix  work  further 
suggested  that  a  major  isomer  of  this  structure  is  a  branched  ring  compound  0  -  Al^? .  This  is  the  same 

configuration  as  an  SCF/CI  calculation  has  indicated  for  the  BO;  intermediate  of  eq.  9.29  aIO  has  been 
observed  to  form  gaseous  addition  complexes  with  other  molecules.  Thus,  chemiluminescense 
observations  in  the  visible  on  the  A1  +  CO2  +  O3  reaction  system  have  given  evidence  for  the  formation 
of  OAICO2.33  Also,  room  temperamre  kinetic  measurements  have  found  that  the  reactions  of  AIO  with 
CO2  and  C7H4  have  pressure-dependent  pseudo-second-order  rate  coefficients. 3^  These  CO2  and 
C2H4  observations  show  an  interesting  difference  with  the  present  MO  +  O2  findings,  in  that  they 
display  a  non-zero  intercept,  indicative  of  a  contribution  of  second-order  process  at  room  temperature. 
That  work  was  done  at  bath  gas  N2O  pressures  of  270  to  800  mbar.  Since  in  the  present  work  at 
pressures  an  order  of  magnitude  lower,  with  a  less  efficient  bath  gas,  Ar,  the  second-order  component 
did  not  become  evident  till  about  1000  K  this  may  suggest  a  rather  stable  OAIO2  adduct. 
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TABLE  II 


Summary  of  Third-Order  Rate  Coefficients  for  the  AlO/02/Ar  Reaction  System. 


T  ±  ot 
(K) 

Na 

P  range 
(mbar) 

[M]  range 
(10^^  cm'^) 

^2,0  —  ^k2,0 
(cm^  molecule"2  s'^) 

305  ±6 

18 

6.8  -  27.9 

1.6  -  6.6 

5.1  ±  1.7(-30)b 

547  ±  19 

14 

7.2  -  32.4 

0.9  -  4.5 

9.1  ±0.5(-31) 

826  ±  24 

15 

11.1  -  102.6 

1.0  -  8.9 

5.7  ±  0.3(-31) 

1010  ±26 

30 

11.3  -  81.7 

0.8  -  5.8 

3.0±0.4(-31) 

a)  N  is  the  number  of  k,  measurements  at  this  temperature. 

b)  Should  be  read  as  (5.1  ±  1.7)  x  10-30  gni6  molecule-2  s'^ 
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Summary  of  Rale  Coeff.c.eni  Measurements  on  the  AlO  +  reaction  taken  above  1 2(X)  K 
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FIGURE  CAPTIONS 

FIG.  1.  Plot  of  kj  vs.  [M]  for  the  AlO/02/Ar  reaction  system  at  T  =  305  K 
( O )  individual  rate  coefficients 
( - )  Troe  fit  of  measurements,  eq.  3 

FIG.  2.  Plot  of  ki  vs.  [M]  for  the  AlO/02/Ar  reaction  system  at  T  =  547  K 
( O )  vaporization  of  A1  and  5%  O2  used  to  produce  AlO 
( + )  vaporization  of  A1  and  0.5%  N2O  used  to  produce  AIO 
( □ )  microwave  discharge  of  TMA  and  0.5%  O2  used  to  produce  AIO 
(  -  )  Proportional  fit  of  measurements,  eq.  4 

FIG.  3.  Plot  of  kj  vs.  [M]  for  the  AlO/02/Ar  reaction  system  at  T  =  826  K 
( O )  Measurements  taken  using  the  C^H-X^S  transition  with  LIF 
( □ )  Measurements  taken  using  the  transition  with  LIF 

(  -  )  Proportional  fit  of  measurements,  eq.  4 

FIG.  4.  Plot  of  kj  vs.  [M]  for  the  AlO/02/Ar  reaction  system  at  T  =  1010  K 
( O)  Measurements  taken  using  O2  (99.98%)  from  Matheson 
( □  )  Measurements  taken  using  O2  (99.6%)  from  Linde 
(  -  )  Linear  fit  of  measurements  eq.  5 

FIG.  5.  Plot  of  the  third-order  rate  coefficient  data  of  the  AIO  +  O2  reaction.  Uncertainties  shown  are 
±2cr  precision. 


(A)  Current  study 
(  -  )  Best  fit  to  current  study,  eq.  6 
(•)  Calculated  from  ref  10 


FIG.  6.  Arrhenius  plot  of  the  second-order  rate  coefficient  data  of  the  AlO  +  O2  reaction. 
( A)  Individual  rate  coefficients 
(-)  Arrhenius  fit  to  data,  eq.  7 
(•)  Recommendation  from  ref  9 
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